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FOREWORD 


This investigation was sponsored by the Office of Materials Processing in 
Space, National Aeronautics and Space Administration, Washington, D.C. Work 
wa9 carried out and is reported in two parts. Part I reports on the investigation of 
the stability of floating liquid columns confined between rotating disks. Part II 
reports on the development of a stereoptical-tracer particle technique for measur- 
ing velocity fields within transparent fluids and its application to these floating 
liquid column systems. 

Dr. A. A. Fowle was the Director and Principal Investigator. Ms. C.A. Wang 
carried out the experiments, made the observations, and reduced and organized 
the data of Part I, Mr. L. Soto provided the same services for the Part II 
investigation and, in addition, was responsible for the apparatus design and test 
protocol used in the measurement of flow fields. Mr. P.F. Strong was the author 
of the theoretical analyses of resonance in axially vibrated cylindrical columns 
and for predicting the location of the tracer particles from stereo photographs. 
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EXPERIMENTAL STUDY OF THE STABILITY AND 
FLOW CHARACTERISTICS OF FLOATING LIQUID 
COLUMNS CONFINED BETWEEN ROTATING DISKS 


ABSTRACT — Part I 


A low Bond number simulation technique was used in conducting a number of experiments 
to establish the stability limits of cylindrical and conical floating liquid columns under conditions 
of isorotation, equal counter-rotation, rotation of one-end only, and parallel-axis offset. Tests 
which establish the conditions for resonance in cylindrical liquid columns perturbed by axial, 
sinusoidal vibration of one end-face are also reported. All tests were carried out under isothermal 
conditions with water and silicone fluids of various viscosities. 

With few exceptions, test data for rotating columns show that the onset of instability is in 
the C-mode and is correlated very well by a theoretical prediction based on isorotation. The onset 
of instability for a cylindrical column of silicone fluid having a viscosity 100 times that of water 
under equal counter-rotation and single end-rotation is in the first axisymmetric mode. An 
explanation for the rather surprisingly simple C-mode correlation is provided by the results of 
flow visualization studies. The data on the resonance of cylindrical columns are correlated very 
well by a prediction based on a newly developed theory. 


I. INTRODUCTION 


The near-zero gravity environment of an orbiting space station offers unique opportunities 
for improving the process of growing crystals in high-temperature, floating molten zones. In such 
an environment, buoyancy-induced convection, which is difficult to predict and control, is greatly 
reduced, and constraints imposed by gravity on the dimensions of a floating zone are relieved. As 
the allowable stable length of a molten zone in a microgravity env ironment is greater than that on 
Earth, the problem of producing a temperature field within the melt that minimizes composi- 
tional variations and residual stresses in the growing crystal is eased. 

Rotating the crystal and feedstock is a common way of effecting in-depth heating of the melt 
zone via forced convection currents and of reducing the azimuthal temperature gradients that 
would exist in its absence. Although rotation may be used to advantage to induce convective heat 
transport to the inner regions of the melt, unwanted temporal and spatial flow variations may 
occur, depending on the speed of rotation, and an upper speed limit is set by its stability. 

As a general rule, then, in growing single crystals by the float-zone method, one may wish to 
employ rotation and the largest possible length-to-diameter ratio of the molten zone. The limits 
of possible length-to-diameter ratios and rotational rates are set by the stability of the liquid 
zone. Resonances due to vibrational disturbances are to be avoided also. It is within these limits 
that geometric and rotational variations are sought to produce crystals approaching perfection. 
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II. SUMMARY 


A review of the literature containing information useful in determining the stability of liquid 
columns having geometries and operating under conditions typical of the classical float-zone, 
crystal-growing method is presented in this report. In addition, the test methods and the results 
of a number of experiments conducted to establish the stability limits of floating liquid columns 
of cylindrical and conical geometry under conditions of isorotation, equal counterrotation, rota- 
tion of one end-face only, and parallel axis offset are given. Investigation of the conditions for 
resonance in cylindrical liquid columns perturbed by axial, sinusoidal vibration of one end-face is 
also reported. 

All experiments were carried out under isothermal conditions and employed liquid columns 
of small dimensions to make the Bond numbers appropriate to each liquid zone of the order one. 
Both distilled water and silicone oil of various viscosities were employed as test fluids. 

The theories that have been advanced to predict the stability of the columns in various 
modes of excitation are presented in Section III, and predictions based on these theories are 
compared with the results of the experiments completed under this contract and with those 
published previously. In addition, a new theory has been developed to correlate the results of tests 
leading to resonance modes of cylindrical liquid columns under normal vibration. 

In all but a few tests employing rotation, it was observed that the columns first encountered 
instability in the C-mode. This result was consistent with theoretical arguments which predicted 
that this was the most unstable mode, that is, it occurred at the lowest rotational speed for a 
given column geometry, but the rather surprising result was that these tests were correlated 
extremely well by the theoretical expressions derived for conical and cylindrical columns for the 
case of isorotation. This was true regardless of the type of rotation, offset, or magnitude of 
viscosity in the range of 1 to 10 centistokes. In some cases, the columns survived the first C-mode 
only to enter and fail in a second C-mode at a higher driving speed. 

Flow visualization studies provided some insight into this result, if not a full understanding 
of it, up to the point of the onset of the first C-mode instability. For the equal counterrotations! 
and single-end rotational cases, the typical Proudman-Taylor cellular structure persisted (modi- 
fied by the appearance of surface vortices for the case of equal counterrotation at sufficiently high 
rotationrl speeds) up to incipient instability. At slightly higher rotational speeds, orbiting of the 
column commenced. The orbiting circumstance was induced when the rotational speed of either 
end -face corresponded to the Rayleigh condition of equality between the change in kinetic energy 
of the column in orbit and the change in potential energy due to surface distortion. This condition 
apparently prevailed irrespective of the circulation within the columns at incipient instability; 
that is, the superposition of fluid force effects held. 

The exceptions to the onset of instability in the C-mode involved tests of cylindrical 
columns of 100-» °ntisoke silicone under conditions of equal counterrotation and rotation of one 
end only. Here the C-mode was not observed and the columns first distorted and failed in their 
first axisymmetric instability mode. It was only in tests involving a viscosity 100 times that of 
water that this effect on stability became noteworthy and it acted to increase the stable range. 
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Agreement between the theory developed to predict the natural frequencies of cylindrical 
columns of water and l-centistoke silicone under normal vibration and test results proved to be 
good. The natural frequencies of columns of silicone having viscosities of 10 and 100 centistokes 
were not evident, the system apparently being overcritically damped or sufficiently damped to 
avoid the observation of resonance. 
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III. THEORY 


General 

The system whose limits of stability are sought is an isotnermal liquid column suspended 
between two rotating disk*. It is characterized by low Bond numbers; that is, the dimensions of 
the liquid columns are made small to minimize the effects of gravity compared to capillary 
effects, thus simulating a microgravity environment. In the simplest configuration, the column is 
a right circular cylinder having a diameter equal to that of the end-disks, and the planes of the 
disks are perpendicular to their axes of rotation which are collinear. 

Upon rotation of the disks, shear forces are generated at the end-faces of a column of viscous 
liquid which, in turn, induce circulation within the column. The surface of the column distorts in 
response to this circulation in an attempt to contain it. In general, the distorted shape of the 
column has a larger area and, hence, larger surface energy. Moreover, a constraint on the 
distorted shape of a stable column is set by continuity; that is, the volume of the distorted column 
must equal that of the undistorted column. The column is contained at rotational rates up to the 
limit where the change in the kinetic energy of the fluid equals the increase in potential energy 
represented by the increase in the product of the surface area and surface tension coefficient. At 
higher rotation rates the column becomes unstable and may fly apart. 

Ix>oked at from the standpoint of a force balance, the surface will distort such that its 
curvature divided into the surface tension coefficient of the liquid will equal the discontinuous 
pressure difference across the surface every place on the surface in order for stability to be 
maintained. In the general case, the flow field within the column is complex, and the curvature of 
the steamlines in both axial and radial planes contributes to the pressure at the surface. However, 
in most cases of interest, the tangential components give rise to the dominant centripetal forces 
which must be balanced by the surface tension forces. As these centripetal forces are proportional 
to the square of the angular velocities of the circulating liquid, a limit on the allowable angular 
velocities of the end-faces exists wherein the distortion of a given column can provide the surface 
tension force required to contain in the liquid and still meet continuity requirements. 

Another form of instability arises when disturbances to the column give rise to capillary 
waves on the surface. The amplitude of these waves is enhanced at resonance modes having 
frequencies occurring when the column length corresponds to an integer number of half-wave 
lengths of the capillary wave. 

Resonances occurring due to normal vibration of the end-faces represent another circum- 
stance of interest. The conditions to be met by resonant modes are that the volume of the fluid in 
the distorted column be constant and that the change in surface energy and the change in kinetic 
energy of the fluid volume be equal. 

The application of dimensional analyses provides insight into the form of the solutions to 
the stability problems under investigation. Equation (1) represents the general form of all 
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solutions under the conditions postulated and where, in addition, viscous forces are insignificant. 
This equation follows from application of standard procedures of dimensional analysis: 


> f (2, 

pajcu, a t u> 2 a 2 


0 ) 


where 

a = surface tension coefficient of the liquid in its surroundings (dyne/cm), 

p - density of the liquid (g/cm 3 ), 

cj x = angular velocity of disk at one end-face (rad/sec), 

cj 2 = angular velocity of disk at other end-face (rad/sec), 

fi = length of liquid column (cm), 

a, = radius of liquid column at one end-face (cm), and 

a 2 « radius of liquid column at the other and-face (cm). 

a 

The term pa 3 co* represents the ratio of surface vision to inertia forces, and is the 

inverse of the Weber number. The physical reasoning cited above provides the basis of the 

inequality sign for stability. The equal sign sets the stability limit. 


In any real system, deviations from the ideal geometry are present, depending on the 
precision of the construction of the assembly. The axes of rotation may be parallel but not 
collinear (offset), or the axes may be non-parallel. In either case, the end-faces may not be 
precisely perpendicular to the rotational axes. In principle, these dimensional variables enter as 
additional dimensionless indices in Equation (1). 


The theories that have been published do not give solutions for stability in the general case 
of a cylindrical liquid column under eccentric or non-parallel conditions, or even to Equation (1), 
when <a)i and differ in magnitude and sign. However, solut ons of simpler cases which provide 
guidance to the design and operational parameters necessary for stable behavior do exist. 
Fortunately, in most cases, practical designs can reduce tho effects of machine tolerance to 
insignificance, as will be demonstrated by test results to be presented. The results of analyses 
pertinent to the major instability problems of practical interest are presented in the following 
paragraphs. 


C-Mode Instability 

Figure la shows the configuration of a distorted liquid column excited in its lowest order 
“C” mode. The C-mode is descriptive of its distorted shape. For steady-state conditions, where 
the end-faces are isorotational, the column rotates as a solid body. In the ideal case of no 
eccentricity and parallel end-faces, the lowest order condition for stability, as developed in 
Reference 1, is given by 

*c = ~J-2 > (-) (2) 

c pa 3 cu 1 \ira/ 
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As before, the equal sign substitutes for the inequality sign at the stability limit. 

The criterion represented by Equation (2) is plotted in Figure 2. Stable conditions of 
op*>«; are those represented by the region above and to the left of the line, <t>c • Unstable 
♦v ui •. (i 're represented in the region below and to the right of this line. 

Equation (2) establishes the maximum rotation rate for a column of given geometry; i.e M a 
right circular cylinder for stable operation under isorotational conditions. In many cases of 
interest the assumed conditions may not be present. The effect of some variations of practical 
interest as determined by experiment will be presented in a following section. 

An approximate expression for the C-mode instability for a conical zone in isorotation is: 
a > [ 0.572a 2 + 0.856a + 0.572 j / g V (3) 

paiw 2 (a+D J \ffa 2 / 

where 

a = a,/a 2) 

a, = radius of simpler end of conical zone, and 

a 2 = radius of larger end of conical zone. 

This expression was derived in a manner similar to that described in Reference 1 for application 
to the case of a cylindrical zone. The maximum value of the bracketed term is 1.0 and occurs for a 
value of aj/a, equal to unity. Therefore, the least stable configuration is a cylinder of radius aj. 
The cone is increasing in stability as u becomes smaller, with the minimum value of the 
bracketed term approaching 0.572 as a approaches zero. Moreover, Equation (3) predicts that the 
stability of the cone is intermediate between a cylinder of radius a, and a 2 , but nearer that of the 
larger cylinder. 


Axisymmetrical Mode Instability 

Figure lb illustrates the distorted shape of the column excited in its lowest order axi- 
symmetrical mode. Hocking 1 derives the stability criteria for this case from an analysis of 
infinitely long liquid columns rotating as a solid body and subjected to axisymmetric longitudinal 
disturbances of small amplitude induced by centripetal forces with the result: 


^xs 


p a 3 to 2 


- 1 


where 


(4) 


2ira _ 2 iran 

X "IT ’ 
a = radius, and 

X = wavelength of disturbance. 
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In the case of columns of finite length, 6 , the instability disturbance be a standing 
wave which does not alter its volume so chat the wavelength must meet the condition *\ - n8 , 
where n is an integer. Therefore, the lowest order condition for stability becomes, for n = 1: 



This criterion establishes the maximum rotation rate for a cylindrical column for stable 
operation under the conditions assumed, and it is illustrated in the curve of Figure 2. Coriell et 
•l.'*’ have analyzed the case of a cylindrical and conical column, including the effects of gravity, 
variant volume, and isorotation. The solution requires machine computation and the results of 
only a few cases have been published. Again, in many cases of interest, the condition of 
isorotation may not be imposed. 


Capillary Wave Resonance 


The generation of longitudinal capillary waves on the surface of a cylindrical column have 
resonant standing wave modes correlated by: 



or, in simp'e~ < -m 


^CW ‘ C 3 W 2 " n 3 ff 3 

Equation (6) represents the simplest version of the many standing wave modes that are possibly 
generated by a variety of disturbances. Here n is the number of half wavelengths in the distorted 
length £ and can take on values 2,3,4, . . . etc., corresponding to s jessiv ely higher modes of 
resonance. A value of n equal to unity is excluded because it violates the condition of continuity of 
volume. Equation (6) with n set equal to 2, corresponding to the first natural frequency, is shown 
as a curve in Figure 2. 

The capillary wave resonances do not involve rotation and consequence centripetal forces as 
a fundamental ingredient, thereby differing from the C-type and axisymmetric modes. Lamb 41 
gives a classical discussion of the physical characteristics of capillary waves. 


Column Resonance undo; Axial Vibration 

The resonant frequencies of axially vibrated cylindrical columns have been predicted on the 
basis of a theoretical analysis developed in the Appendix. At resonance, the free surface is 
distorted by a standing wave containing an integer number of half wavelengths, n, in the 
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length, £ A value of n equal to unity is excluded because it violates constraints set by the 
continuity of the distorted volume. For even values of n, the resonant modes are correlated by 


0av ” p£ 3 (u J = U(P) n 3 l 3 


For odd values of n, the resonances are given by 


(7) 


0av 03 2 ~~ 3 3 

d plror rrir 

The expressions u(p) and v(p) are rather complex functions of p = nm/C 
functions) given in the Appendix. 


( 8 ) 

(involving Bessel 


Summary 

The theor> r relating to the instability of cylindrical columns of liquid is extracted from a 
selection of those that appear in the literature. The scope of the survey has been established by its 
purpose which is to give guidance to the prediction and interpretation of instability phenomena 
occurring in systems of practical significance in crystal -growth processes carried out in a micro- 
gravity environment. The theoretical equations for the first C -mode [Equation (2)], axisymmetric 
mode [Equation (5)), and capillary wave [Equation (6)| instabilities are shown as curves on 
Figure 2. An obvious observation to be drawn from Figure 2 is that the C-rnode instability is 
expected to be the first encountered (the most unstable), and systems designed to avoid this 
instability should, a priori , avoid the others. The theory developed for the natural frequencies of 
an axially vibrated cylindrical zone is believed to he original. In the following sections, experi- 
mental results are presented which test the validity of the stability theories given and their value 
as design criteria. 
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IV. REVIEW OF PRIOR TESTS 


General 

In this section, results of stability tests taken from three sources in the literature are 
reviewed and compared with the theories given previously. These sources are: the work of 
Carruthers and Grasso, 51 the experiments conducted in Skylab IV,*’ and the tests of Fowle et al. tl ‘ 


Experiments Using Plateau Simulation 

Carruthers and Grasso used the Plateau method to investigate the stability of liquid 
columns. Their apparatus consisted of a vertical column of a mixture of water and ethanol 
suspended between vertical, rotatable end-support tubes immersed in mineral oil (Nujol) with 
which it is immiscible. The end-support tubes were hollow to facilitate the deployment of the 
columns. The composition of the vertical column was adjusted in density to equal that of the 
Nujol surroundings as close as practically possible to simulate operation at microgravity condi- 
tions. This apparatus is shown in Figure 3. 

In this work, the investigators studied and reported on column stability under various rota- 
* tional conditions with end-support tubes of equal diameter and different diameters (conical 
zones). The results of tests on the stability of spherical liquid zones, as well as flow v sualization 
studies on cylindrical columns using dye tracers, were alsc reported. 

In the tests on conical and cylindrical columns, isorotational conditions were imposed. The 
axisymmetric mode of instability was observed and correlated with Equation (5) for n in a range 
1 < n ^ 2 for cylindrical zones. Conical zones, formed by using tubes with unequal end- 
diameters, were also tested under co-rotation conditions. Results indicated that the onset of the 
axisymmetric instability occurred as if the conical zone behaved as a cylindrical zone, with the 
limit correlating with the diameter and speed corresponding to the end-face, resulting in the first 
appearing instability. 

The Nujol surroundings suppressed the C-mode instability and prevented solid-body rota- 
tion within the zone, ev n under isorotational conditions, because the boundary conditions on the 
column were different from those of a column in the relatively inviscid surroundings common to 
float-zone crystal operations. The authors also reported an important influence on the tube 
boundarv conditions, whether open or closed, for their Plateau method experiments. 


Skylab IV Experiments 

The apparatus used in the Skylab IV tests is illustrated in Figure 4. Liquid zones of water, 
water plus dye, soap solution, and soap foam were deployed between two aligned circular disks 
which were rotated in isorotation, equal counter-rotation, and single disk rotation only. This 
apparatus made ingenious use of the dedicated time of the astronaut and the limited equipment 
that were available. 
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The experiments were carried out in a microgravity environment in the range of 10 * g, and 
direct visualization and camera observation of a columnar system with a scale in the order of 
centimeters were possible. The system was limited as regards the precision of its geometry, the 
control of rotational rates, the ability to determine the rotational rates of the fluid zone in the 
axisymmetric instability mode, and the time available in one operational sequence to reach 
steady state. 

All the instability inodes discussed in Section II — 0-mode, axisymmetric, and capillary 
wave - were observed in the Sky lab experiments. The axisymmetric modes were observed to he 
dominant under conditions when only one disk was rotated. The generation of these modes took 
place within a transient period, and no equipment was available to measure zone rotational rates. 
The rotational rates were estimated by equations based on the rate ■>! transfer of momentum hy 
shear forces along the axis of the fluid column. These estimates were only approximate. Never- 
theless. the correlation with Equation (4) using the estimated zov* rotational raie was good for 
engineering purposes. 

The C-mode instability was observed to be dominant for ratios of £/ a in the range of ft to 5, 
for both isorotation and equal counterrotation of the end-disks. The C-mode instability was not 
observed at an <2/ a ratio of about 1.5, The observation of the C-mode instability was necessarily 
made under transient conditions only; however, because of its asymmetric nature, the zone 
rotation rate was dearly measurable. 

No marked consistent differences in the instability limits attributable to the mode of 
rotation, no matter whether water-dve, soap solution, or soap foams were used, could be identi- 
fied within the limits of these experiments. The C-mode instability was correlated using an 
analytical expression since determined to be in error. 

Capillary waves were observed to be a dominant feature for an C/a ratio near 1.5 for both 
isorotation and equal counter-rotation conditions. Again, observations were necessarily made 
under transient conditions. Correlation of the rotational rates at which capillary waves were 
observed with a theoretical prediction of first natural frequencies was made. However, the 
theoretical expression used for the comparisons has since been determined to be in error. 


Experiments Using Low Bond Number Simulation 

Fowle et al. 1 used a iiquid column of small geometric scale to simulate stability conditions 
attendant to operation in a low-gravity environment. Figure 5 illustrates schematically the 
experimental apparatus used. The rotation of the end-faces of the column could be controlled and 
measured in magnitude and direction with precision, as well as the rotational rate of an unstable 
fluid zone. The length of the column and the offset of the parallel axes were controlled and 
measured by micrometer adjustments. The apparatus of Figure 5 was also built to accommodate 
flow visualization studies, hence its offset construction. 

This test set-up had the advantages of duplicating the boundary conditions common to a 
true floating zone and the ability to control and measure system conditions with a relatively high 
degree of precision. Its main disadvantage was the fact that gravity forces were not negligible in 
respect to surface tension forces. The influence of gravity forces were tolerable in that they were 
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perpendicular to the centripetal forces, driving both the C-mode and axisymmetric mode of 
instability. The experiments were conducted at Bond numbers, pgafi/a, in a range between 0.61 
and 2.26. 

Distilled water was suspended vertically between rods having parallel axes of rotation and 
end-faces perpendicular to these axes. The column was produced by introducing water from a 
hypodermic syringe in amounts which resulted in a shape approximating a right circular cylinder. 
The shape actually produced could be likened to a pear truncated at top and bottom, adjusted 
such that the cylindrical surface of the fluid column was tangent to the support rod at the top. 
The distortion from a true cylinder was due to gravity forces and, in degree, dependent on the 
length- to-diameter ratio of the column, with the shape approaching a cylinder at smaller £/a 
ratios. 

Tests under isorotational conditions were made. Only the C-mode instability was observed. 
The agreement between test results and the prediction of this instability limit by Equation (2, 
was excellent. 


Summary 

Experimental data from the reviewed literature have been organized in accordance with 
Equation (1) and are displayed in Figure 6. The theoretical equations for the first C-mode 
(Equation (2)), axisymmetric mode [Equation (5)], and capillary wave [Equation (6)] in- 
stabilities are shown as curves on this figure. 

From this display it should be noted that the data from Reference 1 agree very well with the 
theoretical prediction of first occurring instability in the C-mode. The data of Reference 6 
pertaining to isorotation also show tolerable agreement with thi9 prediction, and the differences 
may well be due to the difficulties encountered in obtaining precise data under the circumstances 
of the Skylab tests. The Skylab tests for conditions of single-end rotation similarly show tolerable 
agreement with theory, and the differences may be similarly explained. It is also of interest to 
note here that the C-mode was not the first occurring instability observed in the case of single-end 
rotation. 

The experimental data of Reference 5 do not agree well with the theoretical predictions of 
the first axisymmetric instability mode which was observed. This disparity is believed to be due 
to the presence of shear on the “free” surface of the rotating column imposed by the Nujol 
surroundings used in the Plateau simulation of zero gravity where the actual and theoretically 
assumed boundary condition is nearly shearless. 
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V. EXPERIMENTAL INVESTIGATION 


Description of Experiments 


Summary 

The experimental apparatus and techniques reported n Reference 1 have been extended to 
investigate stability under a wider variety of condit Tn T‘iese include equal counter-rotation of 
the end-faces, rotation of one end-face only, as v eil if rotation under conditions where the 
parallel axes are offset. Water and silicone fluids wd h vis jsities ranging from 1 to 100 centistokes 
were used to determine the effects of viscosity on stability, if any. Cylindrical zones whose end- 
faces were not parallel were investigated for all of the above mentioned conditions of rotation. The 
stability of liquid zones of a conical geometry was investigated also. Finally, using a simple 
apparatus constructed for the purpose, cylindrical floating zones of water and silicone fluid were 
disturbed by the axial vibration of one end-ftce to d termine conditions for resonance. A 
summary of the conditions tested is given in Table 1 while Table 2 summarizes the fluid 
properties used in the reduction of the experimental dau Further details regarding each type of 
experiment follow. 

Rotation of Cylindrical and Conical Zonet; 

The experimental apparatus (previously described tn Section IV) used to determine the 
limits of stability of liquid columns for various conditions of rotation is shown schematically in 
Figure 5. Improvements in the experimental apparatus and technique previously used were 
implemented to enhance the accuracy and reduce tho scatter of measured data. 

In early runs, different degrees of wetting at thi er.d-taces were observed. The wetting was 
caused by machine-made grooves in the surfaces of t \ e Lucite rods. This problem was eliminated 
by lapping the end-faces against each other. Thereafter, liquid columns were produced with well- 
defined contact lines Fixed by the edges of the end-fat ?s cf the Lucite rods. 

Preliminary tests also showed that stability limits were influenc by the volume of liquid 
between the end-faces for a given spacing. Early failure (,it ’ower rotational speeds) occurred for 
liquid zones who r ‘ volume was somewhat greater than that equal to a right circular cylinder, and 
extended stability was observed for liquid zones whose volume was somewhat less than this 
geometry. To minimize the scatter of data due to the uncertainty in the volume of the liquid zone, 
a micrometer syringe was used to inject the volume of liquid equal to that for a right circular 
cylinder with an accuracy estimated to be within 5 r < . Also to minimize the volume changes 
associated with evaporation of the liquid water zones, the z^nc was enclo ed in a Lucite cell which 
contained a reservoir of water at its bottom. Finally, in an attempt to minimize the effect of 
contaminants on the surface tension of the liquids used, the end -/at were cleaned with copious 
amounts of propanol before the liquid zone was formed, 

Cylindrical Zones — Liquid zones were formed with ustilled water and Dow Corning 200 
series silicone fluids. The viscosities ci the silicone fluids were 1, 10, and 100 centistokes. The 
^adius of the cylindrical liquid water zones was 0 4 cm ar.d the length varied to give Bond 

numbers, ggal/a, in the range of 0.68 to 2.05. To achieve comparable Eond numbers for the 
silicone fluid zones (whose surface tension is about one-th.rd that rf water), the radius of the 
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cylindrical zone was reduced to 0.131 cm, giving Bond numbers in the range of 0.92 to 1.8. The 
Lucite rods were coated with an antimigration fluorocarbon resin (Minnesota Mining and Manu- 
facturing's FC-721) to prevent the silicone fluids from creeping up along their sides. 

The effect of offset on stability was investigated for the water zones. The parallel axes of 
rotation were offset by 0.05 cm, corresponding to a 20% radial offset. 

To investigate the effect of non-parallel end-faces, a plexiglass disk, 0.254 cm radius and 
with one inclined face, was glued to the end-face of the botton support rod. The normal to the 
non-parallel end-face was inclined 2.86 degrees to the axis of rotation. Distilled water was used in 
tests with non-parallel end-faces. The length of the column was varied to give Bond numbers in a 
range of 0.86 to 2.05. 

Conical Zones — Conical liquid zones were formed by injecting distilled water in amounts 
equal to the appropriate theoretical volumes between the support rods whose end-faces had 
different diameters. The diameter of the bottom end-face was kept constant, 2a a = 0.507 cm, and 
two diameters of the top end-face were used: 2a t = 0.165 cm and 2a! = 0.254 cm. The bottom-end- 
to-top-end radii ratios were, therefore, 3.07 and 2.0 foi the two cases. Means for offset and non- 
parallel end-faces were also incorporated in the tests on conical zones. 

Vibration 

In addition to investigating the stability criteria for rotating liquid columns, cylindrical 
liquid columns were vibrated axially so that free-surface oscillations and resonant frequencies 
could be observed. The experimental apparatus which was used is shown schematically in Figure 
7. The bottom support rod was positioned along the axis of the sonic vibrator, and the top support 
rod was rigid and coaxial. The support rods were made of machined brass or stainless-steel. 

Distilled water or silicone fluid, 1 centistoke only, was injected between the two end-faces 
through a micrometer syringe in a volume appropriate to a right circular cylinder. To minimize 
volume changes due to evaporation, a water-filled Lucite cell enclosed the liquid water column 
and support rods as before. The axial vibrations were generated and controlled by a Hewlett- 
Packard wide-range oscillator and by a Krohn-Hite amplifier. The zone was observed in both full 
illumination and by stroboscopic lighting through a low-power binocular microscope. 

A range of length-to-radius ratios of about 1 to 4 were investigated. The length of the zone 
was controlled by movement of the top support rod by a screw adjustment. Tests with cylindrical 
liquid zones of different radii were made to determine the dependence of resonant frequencies on 
the radius. For water, end-faces having three different radii were used: 0.666, 0.127, and 0.248 cm. 
For tests with silicone, the radii of the end-faces were 0.652, and 0.131 cm. The full amplitude of 
the vertical displacement of the driven end-face could be controlled and was kept in the range of 
40-70 ^m at all frequencies. The vertical displacement of the vibrating end-face was not altered 
by the present of the liquid zone. 


Results 

Rotation of Cylindrical and Conical Zones 

The experimental data on the stability limits observed for cylindrical and conical liquid 
zones driven by their end-faces in various types of rotation have been reduced to the dimension- 
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less stability parameter, a/pa'u 2 , and have been plotted as a function of C/a in Figures 8 
through 22. The theoretical prediction for the lowest order of instability, the first C-mode, is also 
illustrated in these figures. 

Isurotation of Cylindrical Zones — The experimental results for the isorotation of cylindrical 
liquid zones of water and silicone fluid with and without offset are given in Figures 8, 9, and 10. As 
indicated, in all cases the first occurring instability was the C-mode and the test data are 
correlated very well by the theoretical prediction of Equation (2) (see curve 0 C ). The lower limit of 
stability demonstrated by columns having length-to-radius ratios greater than 2 is probably due 
to the alteration of the shape of the stable column by gravity. The decreased stability of the longer 
columns is a common feature of all tests, regardless of the rotational circumstance. 

On the whole, the first occurrence of instability was unaffected by offset isorotation, even 
though the radial offset was 20 r r. In some cases of offset rotation, however, the column survived 
the 1st C-mode distortion. As the rotation rate was increased beyond that of the first C-mode, the 
column returned to its cylindrical shape, but continued to rotate near to, but slower than, the 
frequency of this first mode. At some still higher speed of the end-faces, a second C-mode 
distortion occurred which resulted in the break-up of the column. The rotation rate of the support 
rods at which the second C-mode occurred was about 1.4 times that of the first C-mode, and the 
zone rotation rate at the second C-mode instability was one-half that of the end-faces. 

The results of tests using silicone fluids having viscosities of 1, 10, and 100 centistokes are 
given in Figure 9. The results show that the first occurring dimensionless stability limits for the 
viscous fluids were only about 10 r f greater than those obtained for the water zones having 
comparable Bond numbers. 

The effect of non-parallel end-faces on the stability of water columns in isorotation is 
illustrated in Figure 10. Non-parallelism, with or without offset, apparently has little influence on 
the occurrence of the first C-mode. However, the column never survived the first C-mode for 
rotation under these circumstances. 

Equal Counter-rotation of Cylindrical Zones — The demonstrated instability of cylindrical 
zones of water and silicone, with and without offset and non-parallel end-faces, when acted upon 
by end-faces in equal counter-rotation, are illustrated in Figures 13, 12, and 13. In all cases, 
except one, the onset of instability is correlated very well by Equation (2) which was derived for 
the case of isorotation. An explanation of this rather surprising result is that the column, when 
rotated at speeds approaching that of the first C-mode, “locks-on” to one end-face and the zone 
rotates at the speed of this end-face, thereby duplicating closely its behavior in isorotation. This 
behavior was revealed by flow visualization studies which are described in more detail in a 
following section. 

It should be noted that the stability of the column with counter-rotation of the end-faces is 
very similar to that with isorotation, including the relatively minor influences of offset and 
nonparallelism of the end-faces, although the ability of the column to survive the first C-mode 
was somewhat enhanced. Viscosity in the range of 1 to 10 centistokes has a minor influence on 
results, but at 100 centistokes, no C-mode distortion was observed The single point at 100 
centistokes, shown in Figure 12, represented distortion and failure in the first axisymmetric 
(bottle-shaped) mode. The theoretical prediction for this mode is shown as a curve, 0„, in Figure 
12 . 
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Single-end Rotation of Cylindrical Zones — Test results on the stability of cylindrical zones 
of water and silicone, with and without offset and non-parallel end-faces, when driven by the 
rotation of one end face only, is illustrated in Figures 14, 15, and 16. Again, in all cases except 
those involving the use of 100-centistoke silicone, the onset of instability is in the C-mode and 
results are correlated very well by Equation (2). 

The stability of the columns with single-end rotation proved to be very similar to those 
driven in isorotation or counter-rotation. The reason for this result was revealed by flow vis- 
ualization. The column was shown to rotate pretty much as a solid body at the speed of the 
rotating end-face. Therefore, its first mode of instability was very nearly identical to that of a 
zone driven in isorotation. The relatively minor effects of viscosity (in the range of 1 to 10 
centistokes), offset, and non-parallelism were again demonstrated. The somewhat enhanced 
ability of the columns to survive the first C-mode was similar to its behavior when driven by equal 
counter-rotating end-faces. The observation of the absence of the C-mode and of the oneset of 
instability and failure in the first axisymmetric mode for tests made with 100-centistoke silicone 
fluid was also similar. 

This result is in variance with that observed in the Skylab tests* where the onset of 
instability under similar conditions was in the first axisymmetric mode (see Figure 6). The reason 
for this significant difference is not clear. However, it is noted that the Skylab tests were not 
carried out in steady state, and it may be that the system was rapidly accelerated through the 1st 
C-mode without failure which these tests show to be entirely possible. 

Isorotation and Corotation of Conical Zones — The results of tests on conical water 
columns with isorotation and corotation of the end-faces are shown in Figures 17 and 18. With 
isorotation, the column failed at the first C-mode at a speed corresponding to the first critical 
based on the large end-face (bottom end). The results are correlated better by Equation (2) than 
by Equation (3); that is, the conical zones behaved as if they were a cylinder having the diameter 
of the larger end. The increased degree of stability predicted by Equation (3) did not manifest 
itself. The results of a few tests in which the rotating axes were offset do not alter this picture. 

In corotation with a larger end-face rotating at a speed well below that corresponding to the 
first critical speed, the column was observed to distort in the C-mode when the speed of the 
smaller end-face (top end) reached the first critical speed corresponding to the larger end, but it 
did not fail. The column returned to its normal 9hape with increasing top-end speed up to the 
limit of the rotary speed of the motor, even though this speed might have been greater than the 
critical speed corresponding to a cylinder having the diameter of the top-end. 

In corotation, with the larger e* a rotat'ng at a speed relatively close (say, within 75 r f ) to the 
first critical corresponding to this end, the f r*t onset Oi stability could occur when the speed of 
the smaller end exceeded this critical. The zone rotation under these circumstances had a 
frequency intermediate between the speeds of the two ends above, but not far from, the first 
critical observed in the case of isorotation. 

Equal and Unequal Counter-rotation oj Conical Zones — The results of tests to demonstrate 
the stability of conical liquid zones when driven by end-faces in equal and unequal counter- 
rotation are shown in Figures 19 and 20. Water was uaed in these tests. The results were 
remarkably similar to those of isorotation and corotation, which suggests that the zone rotated 
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much as a solid body attached to the larger end when the speed of the larger end approached its 
first critical C-mode. Again, results correlate well with Equation (2), 


Single-end Rotation of Conical Zones — Results of tests on the stability of conical zones of 
water with rotation of the larger (bottom) end-face and smaller (top) end-face are illustrated in 
Figures 21 and 22, respectively. With larger end rotation, the column behaves as it did in 
isorotation. When the column was driven by the smaller end, its stability limit was set by the first 
C-mode based on the larger end. Also, the column was a little more stable when driven by the 
smaller end. The onset of instability was a bit delayed and the column distorted, but did not fail 
as the speed was increased through this first critical. However, these differences were not very 
significant and the results are correlated well by Equation v 2). 

Axial Vibration of Cylindrical Zones 

Resonant frequencies of axially vibrated liquid columns were observed as standing waves on 
the free surface. Observation under stroboscopic lighting revealed that there were several modes 
of resonance of increasing natural frequency. At resonance, the length of the free surface distorted 
to form a standing wave having an integer number of half wavelengths, with the exception of 1.0, 
which circumstance cannot meet the continuity of volume requirement. The shapes of the liquid 
zone for the first two modes are shown in the photographs of Figure 23. 

The results of all data taken with water and a silicone fluid having a viscosity of 1 centistoke 
are shown plotted in Figures 24 and 25, respectively. The resonances were not obvious when 
silicone having viscosities of 10 and 100 centistokes were used. The data are plotted in accordance 
with the correlations suggested by the theories developed in the Appendix, and curves represent- 
ing the theoretical predictions are superimposed. 

The comparison between experiment and theory is, on the whole, useful if not precise. The 
differences, particularly in the case of the odd number modes, may well be attributable, in part, 
to the appreciable amplitude of the perturbation that was applied where the theory assumed an 
infinitesimal disturbance. The disparity may also be, in part, due to the influence of gravity, 
which plays a larger role in the cases involving the larger ?/ a ratios. Nevertheless, the differences 
are not large when viewed from an engineering perspective, and the forms of the functions u(p) 
and v(p) are in good agreement with the observations. 


Visual Observations 


Method 

As previously stated, the onset of instability was determined by visual observation of the 
shape of the liquid zones. These shapes were made evident by illumination from either a high- 
intensity, tungsten-halogen, incandescent lamp or a stroboscopic, xenon-flash light source. Fig- 
ure 26a shows the distortion of the zone in the first C-mode and Figure 26b shows its shape just 
prior to failure. 

A limited number of tests were also made to identify the circulation patterns within the 
liquid zones as they were rotated in the approach to their stability limits in order to provide better 
insight into the physical reasons for the observed results. For this purpose, nickel plated, hollow 
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glass spheres, 5 to 50 nm in diameter, were introduced into the column. These spheres were 
previously screened for neutral buoyancy by letting them settle in a column of water. All flow 
visualization tests were made on a cylindrical water column 0.254 cm in radius and 0.510 cm high. 
The motion of the tracer particles were observed with scattered light from a dark-field illumina- 
tion system. The flow patterns in the liquid columns were observed under conditions of isorota- 
tion, counter-rotation, and single-end rotation. 

Results 

When the column was driven in isorotation it behaved like a solid body as expected. The 
tracer particles were observed to whirl in rings with no evident radial or axial velocity component. 
When the water column was rotated in equal counter-rotation at frequencies low in respect to the 
first C-mode critical, two Proudman-Taylor cells were generated, one on each side of an inter- 
mediate axial plane. As the rotational rate was increased, the bottom cell became larger, no doubt 
because of the influence of gravity on the shape of the zone. At all rotational speeds, the flow 
spiralled toward the disks in a central core, then spiralled radially outward at the end-faces, 
thence spiralled away from the disks in an annular layer next to the free surface, and completed 
the circuit by spiralling radially inward at an intermediate layer. The reversal of the tangential 
velocity component was accommodated in a relatively thin Ekman layer in an intermediate plane 
between the end-faces. 

Another interesting feature of the circulating flow wus also seen. At about 200 rpm, two or 
three vortex swirls, with radial vortex lines, appeared on the free surface at the joining region of 
the two Proudman-Taylor cells. These have been reported previously in Reference 1. These 
vortices moved upward to the top end as the speed was increased. Their centers moved in the 
direction of the top end-face with some small relative tangential velocity. 

With single-end rotation approaching the critical speed, the bulk fluid circulated with a 
large Proudman-Taylor cell associated with the rotary speed of the driving end-face. Again the 
flow spiralled down the central core away from the rotating end and returned in a spiral path in 
an annular zone next to the free surface. No surface swirls were seen for this type of rotational 
drive. The radial inflow and the reversal of the tangential component of velocity took place in an 
Ekman layer next to the stationary end-face. 

In all cases of single-end or counterrotation, the orbiting circumstance of the first C-mode 
was induced when the rotational speed of either end-face corresponded to the Rayleigh condition 
of equality between the change in kinetic energy of the column in orbit and the change in 
potential energy due to surface distortion. 
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TABLE 2 


FLUID PROPERTIES 


Vitoo«ity 



Dantity 

Kinematic 

Dynamic 

Fluid 

(gm/cc) 

(c.».) 

(dyna-Mc/cm 3 ) 

Water 

1.00 

1 

0.010 

Silicone Fluid 

0.82 

1 

0.0082 


0.94 

10 

0.094 


0.07 

100 

0.970 
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Surface 

Tension 

(dyne/cm) 

72.6 

17.4 

20.1 

20.9 
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Figure la 


C-Mode Deformation 


Figure 1b 


Axisymmetric Deformation 


FIGURE 1 


UNSTABLE SHAPES OF ROTATING LIQUID ZONES 
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Zone Lengthto-Radius Ratio, (/a 


FIGURE 2 


THEORETICAL STABILITY CRITERIA FOR CYLINDRICAL COLUMN 
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FIGURE 3 EXPERIMENTAL APPARATUS USED FOR FLOATING LIQUID ZONE 
STUDIES BY PLATEAU SIMULATION 
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CIRCULAR 
LAUNCH 
RESTRAINTS (4) 


FLOATING 

ZONE 


OAC SLIDE (4),^ 



I 8" EXTENSION 

4* EXTENSION r 13' PINCH BAR 


EXTENSION 


(b) 


FIGURE 4 EXPERIMENTAL APPARATUS FOR LIQUID FLOATING ZONE 
STUDIES UN SKYLAB IV 
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Uppar Driva 
Motor-Controllar 



FIGURE 6 


Illuminator 



EXPERIMENTAL APPARATUS USED FOR FLOATING ZONE 
STUDIES BY LOW BOND NUMBER SIMULATION 
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Stability Limit Parameter 
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Screw Adjustment 



FIGURE 7 SCHEMATIC OF EXPERIMENTAL APPARATUS USED FOR 
AXIAL VIBRATION OF CYLINDRICAL LIQUID COLUMNS 


27 


Arthur D Little Inc 





FIGURE 8 CYLINDRICAL WATER COLUMN - ISOROTATION 
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Stability Limit Parameter, 



FIGURE 10 CYLINDRICAL WATER COLUMN - NON-PARALLEL END-FACES. ISOROTATION 
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Stability Limit Pr 



FIGURE 11 CYCLINDRICAL WATER COLUMN - EQUAL COUNTER ROTATION 
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FIGURE 12 CYLINDRICAL SILICONE COLUMN - EQUAL COUNTER-ROTATION 
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Stability Limit Parameter 



FIGURE 13 CYLINDRICAL WATER COLUMN - NON-PARALLEL END-FACES, EQUAL COUNTERROTATION 
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Stability Limit P> 
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FIGURE IB CYLINDRICAL SILICONE COLUMN - SINGLE END ROTATION 
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Stability Limit Parameter, 



FIGURE 16 CYLINDRICAL WATER COLUMN - NON PARALLEL END-FACES, SINGI E-END ROTATION 
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Stability Limit Parameter, 



0.1 1.0 10 

Zone Length-to-Large End Radius Ratio, 8/«2 


FIGURE 17 CONICAL WATER COLUMN - ISOROTA'ION AND CO-ROTATION 
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Stability Limit Parameter. 



FIGURE 18 CONICAL WATER COLUMN - ISOROTATION AND CO-ROTATION 
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FIGURE 19 CONICAL WATER COLUMN - EQUAL AND UNEQUAL COUNTER-ROTATION 
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Stability Limit Parameter, 



FIGURE 20 CONICAL WATER COLUMN - EQUAL AND UNEQUAL COUNTER-ROTATION 
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Stability Limit Parameter, 





FIGURE 21 CONICAL WATER COLUMN - BOTTON END ROTATION 
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Stability Limit Parameter, 



FIGURE 22 CONICAL WATER COLUMN - TOP-END ROTATION 
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First Mode (n = 2) 




Second Mode (n = 3) 

FIGURE 23 SHAPES OF CYLINDRICAL ZONES AXIALLY VIBRATED IN FIRST 
TWO MODES OF RESONANCE 

I 


43 


Arthur Dl-ittlelrx 





Rtdiu* of Cylindrical Zona 

• a *0.066 cm 
o a -0.127 cm 
■ a *0.249 cm 








Eq. (A— 41) 


Rayleigh Static 
Stability Limit 




Dimensionless Resonance Parameter, 



FIGURE 26 AXIALLY-VIBRATED CYLINDRICAL COLUMNS OF SILICONE 
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di Liquid Zone Rotates in C Mode Deformation 


b) Ltquid Zone Failure in C Mode 

FIGURE 26 CYLINDRICAL WATER ZONE IN C-MODE DEFORMATION 
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APPENDIX 


THE PREDICTION OF THE NORMAL VIBRATIONS 
OF A CYLINDER OF LIQUID 


The Normal Vibrations c e a Cylinder of Fluid 

In Figure A-l we show a circular cylinder of fluid suspended by capillary force' 
between two rigid circular plates. The fluid represents the molten zone in a crystal- 
growing process, and the circular plates represent the ends of the feed rod and 
growing crystal. We assume that the effects of gravity can be ignored, either because 
of the small size of the fluid zone or because the experiment is performed in a weight- 
less environment. 


In Figure A-2 we show a mode (characterized by mode number n=l ) which is not 
possible to realize in the case of a cylinder of fluid because the volume is not conserved. 
Because this mode would be allowed in the case, say, of an elastic membrane stretched 
between two rings to form a hollow cylinder, we give it mode number n=l , although it 
is suppressed. As will be seen later, higher odd modes are present. 


In Figure A-3, we show the mode n=2. The perturbed radius is taken to be 

2jtz 

r = a + a sin — — 

£ 


(A-l) 


where a is small compared to a. The volume is 

V = / irr 3 dz = xa 2 £ + 2iraa j'\ 


= J irr 2 dz = xa 2 £ + 2iraa f sin — dz + 0 (a 2 ) 


(A-2) 


Since the integral is zero, the volume is conserved at least to first order. For the 
even modes, in general, we may write the perturbed radius as 


r 


a + asin 


nnz 

£ 


, n = 2, 4, 6, etc. 


(A-3) 


Next we calculate the potential energy, P, due to the increase of surface area 
S and the surface tension a. The surface area is given by the integral 




2irrds 


(A-4) 
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where ds is an element of arc as shown in Figure A-3. Now 

ds = (dz l + dr 2 )* = £l + % dz 


so, approximately, 

ds '[ l+ *(‘^') ] dz 


(A-5) 


(A-6) 


Therefore the surface area is seen to be 




n 2 7r 2 u: 2 

2£ 2 


COS ' 1 


mrz 


dz 


(A-7) 


Of the four terms in Equation (A-7) obtained by multiplying the binomials, only two 
have averages different from zero. Thus the surface area of the distorted shape is seen 
to be 


S = 2ira£ + 


n 2 rr 3 aa 2 

2 £ 


(A-8) 


Evidently the first term is the unperturbed area, so we let that correspond to 
the zero of potential energy and find the change about this reference 

_ on 2 jr 3 a , 


(A-9) 


The next phase of the problem is to calculate the kinetic energy associated 
with the vibration, and to do this we have to find the radial and axial velocities, 
v f and v z . To do this, we introduce a velocity potential in the form: 


„ n?rz . /nirr\ 

<p = /?sin — I 0 ^ ^ J (A- 10) 

where I 0 is the Bessel function of imaginary argument. 
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This leads to velocities 


V _ 

' ~ dr 

nrr 

e 

. nirz 

- 0 sin — 1, 

K 


(A-i n 

dsf 

n» 

» "I 1 , 

( n * r \ 


Vj ~ ~dz 

"T 

■ & COS - I 0 

(—) 

(A-l 2) 

From Equation (A-l 1) we see that 

the radial velocity is zero at the top 



and bottom (z = 0 and z-5) but from Equation (A-l 2) we see that v is not zero 
at z=0 and z={. Therefore, the vertical velocity is not consistent with the 
appropriate boundary conditions and is incorrect. Sinc^ the cosine is unity 
at z=0 and z=C. we try the velocity: 

v njr „ / ntrz \ / mrr \ 

•’T'rT'V-lT-) 

in place ot Equation (A-l 2), so the boundary conditions are met. The only 
condition wo lu\e to satisfy . then, is continuity. That is: 


Since 


and 


div v - 0 



(A-l 4) 
(A-l 5) 


(A-l 6) 


(which is Bessel's equation), we see that Equation (A-l 4) is satisfied by the 
velocity as determined by Equations (A-l 1) and (A-l 3). In fact, the velocity can 
be derived from the stream function 


with 

and 


* = 01 , 



ntrz 


-) 


v, = 


dr 


= — ~~ (r*) 


(A-l 7) 
(A-l 8) 
(A-l 9) 
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Next we calculate the velocity at the free surface in two ways to identify 
the quantity 0. To begin with, differentiate Equation (A-3) with respect to time: 


v r (a) = a sin ■ 


njrz 


(A-20) 


and next substitute r=a in Equation (A-l 1 >: 

mr0 / nira 


Therefore we see that 


• nir<3 
a = — - — I< 


or 


0 = ' 


nir I, 


/nira \ 

W 


— > sin 

(A-21) 

8 J 8 

/ nira \ 

V* ) 

(A-22) 

— of 

(A-23) 


The kinetic energy is given by the integral 

K = J'vi p (v t 2 + v 2 ) dV 


taken over the volume of the fluid; in Equation (A-24) p is the density and 
dV is an element of volume. Substitution of Equations ( A-l 1 ) and (A-i 3) into 
Equation (A-24) leads to the expression 




nirz 


, nirz nirz 

— ; 2 cos + 1 


2ir ' dr dz 


The integration with respect to z is quit? simple and leads to the expression 

2C -b |_ • \ K / \ 2 /j 


K = 


(A-24) 


(A-25) 


(A-26) 
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Although the functions in the integrand above are less familiar than sines and 
cosines, the integrals are easily evaluated from the formulas 

/xl’o(x)dx= -f HS(x)-I](x)l +C 

/ y 2 

xl](x)dx = xl 0 (x)l,(x) ip Uj(x)-I?(x)l + C • 


and 


These formulas and a bit of tedious manipulation show that 

n ir 2 |3 2 p a 


K = 


{ I 0 (P)MP) + P[Io(P) - Ii(P)l| 


where we have written 


mra 



for short. After substituting the value of 0 from Equation (A-23) into the 
expression for kinetic energy, we find: 


K 


P a E 2 r 1 q(p) t / Iq(p) 
2n Mp) P \ I* (P) 



In orde* - to obtain the resonance frequencies we set 


a~a e iaJ1 


and appeal to the conservation of energy : 

P + K = const. 


whence, from Equations (A-9), (A-31) and (A-32): 


~ a n 2 ir 3 
28 


to 2 p a 
2n 


8 2 


/ mra\l^ 

U \T/J a 


a 1 e 2,wt = const 
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(A-34) 
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where 


u(p) = 


Iq(P) 

I,(P) 


+ 




(A-35) 


Equation (A-34) can be satisfied only if the coefficient of the time-dependent term 
is zero, so we obtain the result 


n 3 ir 3 n 



(A-36) 


The function u(p) can be evaluated for moderate values of p by reference to tables 
of Bessel functions. Beyond that point the function can be evaluate^ -oproximately 
by means of asymptotic expansions of I 0 and I, : 


Io(P) ^ 


eP 


H 


ypirp \ 8p 128p 2 


I,(P)^ 


yj^P 


(- 


15 


8p 128p 2 


By long division we find 


4SBL„ lt JL. » 

I,(p) 2p 8p 2 


So we find 


u(p) ^2 + 


2p 


(A-37) 

(A-38) 

(A-39) 

(A-40) 


At p=iO, for example, the exact value of u is 2.1 67 • • • , while the asymptotic 
expression gives 2.150; the error is about 0.8%, i.e., it is in the order of 1/p 2 . Some 
useful values of u(p) are given in Table A-l. 
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Another way to put Equation (A-36) which corresponds to the manner in which 
the experimental results are organized is: 



The analysis of the odd modes is not so easy as the analysis of the even modes 
because the mode shape cannot be simply sinusoidal or the volume would not be 
constant throughout the motion. 

At this juncture it is best to redefine the z-coordinate so that the origin is at 
the center of the fluid zone, which now is limited at z = ± C/2. As a result of this 
change, the roles of sines and cosines are reversed, but the additional symmetry 
makes the algebra significantly simpler. 

In this case we may start with the stream function: 

i> = 01,(^) (/a sin + a 3 z 3 + a, z) . (A-42) 

Here the coefficients a 3 and a, are to be determined to satisfy the boundary 
conditions and conserve mass, while the coefficient /i is defined as follows: 

M = +1 for n = 3, 7, 1 1. ... 

H = -1 for n = 5, 9, 13, ... 

This is done so that n sin - 1 for all odd values of n; otherwise the two sets of 

odd numbers would have to be treated separately. 

As before we take 

9 6 nnp , /n irr\ / nirz 

v ' " ” 5z r Ht) — 


3 (a. 


nir 


L 2 + 


nir 


(A-45) 


(A-43) 
( A-44) 
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The boundary condition to be imposed is v r ^z = ± = 0. Since cos ~ - 0, 


we 


must have 


^-a 3 + =0 , 

4mr n7r 


or 


a, = 


4 


which result can be substituted into the stream function to yield 


Next we consider v z : 


1 8 nnp . /nirr\/ mrz f, 3C 2 z"]\ 

■7¥ W, *T + 31 [ z -~J) 


Again the boundary condition is such that v 


1 + 


■M 


= 0. So we require that 


or 


a 3 = - 


C 3 


Finally, therefore, the stream function and components of velocity are found to be 


* = JH, 



/i sin 


mrz 


4z 3 
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(A-49) 


(A-50) 


(A-51) 
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d\[/ 

ror/3 

(™\( ucos 

nirz 

12z J 

, + -1\ 

Vr = azT 

T 1 

>vta m 

£ 

nir£ J 

nir/ 

l d 

. ,, nvB 

/n*r\ ( 

nirz 

4z 3 

3z \ 

V ‘ ’ r 3r 

W) = — I 0 

vr)v sm 

~T 

" IP 

+ TV 

The mode shape must have the 

same dependence on z as v r 

Therefore we 


write the mode shape as 

/ nirz 1 2z 2 3 \ 

r = a + a I pcos + — ) 

V £ nirfi 1 nir / 


An illustration of the mode shape for n=3 is given in Figure A-4. The volume is 
then given by the integral 



so the volume is conserved, again to first order. 


To compute the potential energy we start by calculating the change in surface 
area AS from the formula 


AS 



2irr 



dz - 2ira£ 
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The term that is linear in o involves the same integral as in Equation (A-57) and is 
therefore zero. We see, then, that 



which leads to the result that the potential energy P is: 

_ n 2 7r 3 aa 2 o /, 96 \ 

2C ~ V ~ n 5 * 5 ) 

which is nearly the same as it was in the case of the even modes — see Equation (A-9). 

The task of calculating the kinetic energy K is more formidable. We start from 
the basic definition 

/•a * 8/2 , 

K =| I — p (v r 2 + v z 2 )- 2?rr drdz • 

Jo i-m ~ 

On account of the symmetry, we may write: 

r a A/2 

K = 2 7rp I I (v t 2 + v z 2 ) r dr dz 

Jo Jo 

whence 
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(A-62) 


(A-63) 


(A-64) 
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We perform the integrations with respect to z first and discovery, by much effort, 
that 



It is instructive to compare this with Equation (A-26), which gives the kinetic energy ■ 
for the even modes. In that case the coefficients of Io and I? are 3 and I , respectively, 
whereas here they are approximately 2 and 1 . The integrals of the Bessel functions 
can now be carried out just as in the previous case so we obtain the result: 


K = 


nir 2 0 2 pa 


C| i 0 (p)i, (p) 


C 2 p 


[lo(P) - I? (P)]J 


where p = nna/H as before and the constants C, and C 2 (both nearly unity) are: 

„ , 48 192 

C | s l + t~t~i 4 4~ 

Sn 2 * 1 n 4 tr 

34 48 

Cz 35 Jn 2 n r 

The expression for 0 2 is unchanged - see Equation (A*23) — and consequently 


K = 


pafi 1 

2n 


lo(p) Cjp / lg(p) 

- I,(P) 2 \I?(p) 



.2 

a 


(A-66) 


(A-67) 

(A-68) 


(A-69) 


Finally, from the principle of conservation of energy, we obtain the resonance 
frequencies as 


n 3 w 3 o 


(A-70) 
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where v(p) = 


Io(P) 

I«(P) 


C t p 

2 


_ il 

L i?(p) J 


1 - 


96 


I 0 (P) 

Ii(p) 


+ -P 

2 


1q(p) 

1?(P) 


- 1 


This result may be expressed in the alternative form 

a _ 1 /n?ra\ 

pfi 3 w J n 3 ir 3 ' C / 

Table A-2 gives values of the function v(p) calculated from the approximate 
form of Equation (A-71). 
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TABLE A-1 


P 

ir/2 

2ir/3 

it 

4tt/3 

27T 

4ir 

87T 


The function u(p) 


u(p) 

4.239 
3.432 
2.751 
2.489 
2.287 
2.130 
~ 2.060 

|q(p) + / Iq(p) 
•i(p) P \ if(p) 


•) 


TABLE A-2 


P v(p) 

s/2 2.898 

2jt/3 2.396 

it 1.97S 

2 it 1.689 

4 n 1.586 

8it 1.541 

T . . . . UP) p / i3(p) 

The function v(p) ^ + — ! 


li(p) 2 \ l?(p) 


•) 
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FIGURE A— 1 CYLINDER OF FLUID SUSPENDED BETWEEN 
RIGID CIRCULAR PLATES 
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FIGURE A— 2 



THE SHAPE OF THE MODE n - 1 DOES NOT CONSERVE 
VOLUME 
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FIGURE A— 3 THE SHAPE OF THE MODE CORRESPONDING TO n • 2 
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EXPERIMENTAL STUDY OF THE STABILITY AND 
FLOW CHARACTERISTICS OF FLOATING LIQUID 
COLUMNS CONFINED BETWEEN ROTATING DISKS 


ABSTRACT — PART II 


A technique for the quantitative measurement of stream velocity within a floating, isother- 
mal, liquid column confined between rotatable disks has been developed. In the measurement, 
small, light-scattering particles were used as streamline markers in common arrangement, but 
the capability of the measurement was extended by use of a stereopair photography system to 
provide quantitative data. Results of velocity measurements made under a few selected condi- 
tions, which established the precision and accuracy of the technique, are reported. The general 
qualitative features of the isothermal flow patterns under various conditions of end-face rotation 
resulting from both still photography and motion pictures are reported. Also the outstanding 
qualitative features of the flow patterns for liquid columns at low Bond numbers under counter- 
rotational conditions, which may be of special interest to the development of a scientific under- 
standing of the float-zone crystal-growth process, are examined in more detail. 


pRcn: r r"T, t not filmkh 
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I. INTRODUCTION 


An accurate quantitative description of the mass, momentum, and heat transfer within the 
molten zone of a floating-zone, crystal-growth process is necessary to the full scientific under- 
standing of this process. This understanding is a prerequisite for advancing the present tech 
nology which is based on a patchwork of empiricism and partial theory. In recent theoretical and 
experimental work, mainly under the sponsorship of United States and European space pro- 
grams, this need has been recognized, and these programs have added greatly to the state of 
knowledge. However, there is now almost no quantitative te9t data on the behavior of floating 
zone systems of interest to the manufacturer of crystal materials, and a check between developing 
theory and experiment has yet to be accomplished. 

The purpose of the work described in Part II of this study was to provide: 

1. Experimental observations of the flow within floating zones of value to the quest for 
a better scientific understanding of the float zone process, and 

2. A method for measuring velocity within floating zone configurations which can be 
used to obtain a data base against which developing predictive methods can be 
tested and improved. 

The work simulates approximately the low-gravity condition by employing models of small 
geometric scale. 
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II. REVIEW OF PRIOR FLOW VISUALIZATION STUDIES 


A number of investigators have carried out experiments especially devised to observe the 
flow within floating zones under conditions of interest to the materials processing community. 
Carruthers and Grasso' 1 ' made qualitative observations of rotating liquid columns using dye 
markers in a plateau-type simulation of low-gravity conditions, as described in Part I. As noted, 
this method of simulation gave rise to a spurious viscous shear condition at the free surface of the 
column which significantly altered the observed circulation within the column from that which 
would be present with a typical liquid/gas boundary. Furthermore, the patterns observed were 
subject to distortion caused by refraction at the interface of the water-ethanol mixture and nujol, 
as well as the walls of the tank which made their interpretation imprecise. 

Chun and Wuest <JU#Ml have reported a number of useful experiments in which they observed 
the flow patterns within small liquid columns of water, alcohol, and silicone oil with heating and 
end-face rotation using optical/particle tracer and low Bond number simulation techniques. They 
used an optical system which illuminated the meridian plane with a sheet of light. With this 
method, they obtained an uncomplicated picture of the axial and radial velocity field universal to 
the axi-symmetric case of rotation, but no information about the tangential components of 
velocity and inadequate information about possible velocity fields that do not have axial symme- 
try. 


Schwabe et have reported somewhat similar experiments in which they observed 

the flow patterns within small liquid columns of NaNO a made visible by the paths of bubbles 
induced at sufficiently low pressure. Only qualitative visual observations of the flow patterns 
were reported. 

As a precursor to this work, Fowle et al . 8 9) also used low Bond number simulations to 
observe the behavior of floating liquid columns, using the apparatus described briefly in Part I. 
In this earlier work, dark field illumination of small, polystyrene tracer particles was used 
successfully for qualitative flow observations, but the low reflectivity of these particles limited 
their utility to observation of low particle velocities. Also the basic concept for making quan- 
titative velocity measurements using a pulsed flash source and stereopair photography and the 
refinements necessary to develop it to a workable technique were defined. 
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DESCRIPTION OF EXPERIMENTS 


Flow Survey 

The experimental apparatus and techniques reported in Part I were used to survey the more 
general and obvious characteristics of the flow within a floating liquid column induced by 
counterrotating end-faces. One of the flow features which was extensively studied was the effect 
of the velocity ratio of the top and bottom end-faces on the position of the surface of separation 
between Proudman-Taylor cells associated with these end-faces, because this relatively obvious 
feature could provide a simple test of developing predictive analyses. 

Tests were performed with distilled water and Dow Corning 200 series silicone fluids. The 
viscosities of the silicone fluids were 3, 5 and 10 centistokes. Nickel-coated, neutrally buoyant, 
hollow-glass spheres of 35-^m radius were used as flow-tracking particles. A binocular microscope 
with a calibrated eyepiece was used for observation and geometrical measurement. 

The procedure used in making these .measurements was to maintain the angular velocity of 
one of the end-faces constant, while increasing the angular velocity of the other end-face above 
that held constant. The minimum value of the velocity ratio for which data were obtained was the 
value for which the position of the plane of separation became indistinguishediy near to the 
slower end-face, at which point the C-mode instability first occurred, or when limited by the 
rotational speed of the motors. 

As reported in Part I and in previous work, 8 surface swirls occur within cylindrical columns 
having end-faces that are counterrotating at sufficiently high rates. These vortices appear at the 
plane of separation between the Proudman-Taylor cells associated with the upper and lower end- 
faces and have vortex lines radially directed at the surface. Further observations of the flow 
circumstances resulting in these swirls were made, using water as a test fluid for columns having 
radii of 0.254 and 0.127 cm. Both nickel-cue ed, hollow glass spheres and lycopodium powder 
(Sargent Welch item no. SC13245-102) were used to decoiate the swirls. 

The position of the surface of separation of the Proudman-Taylor cells was identified with 
the location of a zero tangential component of velocity at the surface, as revealed by particle 
motion on the surface of the liquid column. For rotational conditions sufficiently high for surface 
swirls to occur, the position of the center of the surface swirls was used. 

It is important to note that the identification of the plane of separation of the Proudman- 
Taylor cells, with a condition of zero tangential velocity at the surface, is not necessarily rigorous, 
This condition was chosen because it was visually obvious and apparently corresponded to the 
surface separating the cells, as revealed by the division between upward- and downward-directed 
axial velocity components. However, because of the uncertainties associated with the observa- 
tion, it was not absolutely clear that the condition of zero tangential velocity at the surface 
corresponded identically with the condition of zero axial velocity, or even that the condition of 
zero axial velocity was precisely at a plane perpendicular to the rotational axis. This issue can be 
resolved by more tests using the quantitative technique for \elocity field mapping described next. 

The conditions tested are summarized in Table 1, anc the fluid properties used in the 
reduction of the experimental data are summarized in Table 2. 
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Velocity Measurements 


The primary objective of the velocity-measurement work was to improve the existing 
technique for flow visualization using tracer particles, to enable the quantitative measurement of 
position and velocity within rotating, floating liquid columns. To accomplish this, the following 
five tasks were undertaken: 

1. The improvement of the light-scattering characteristic of the flow tracer particles. 

This was achieved by substituting nickel-plated, hollow glass spheres for the 
polystyrene spheres used in previous work. These particles were screened both for 
size and density. The particles were sieved to get a cut between 30- and 45-pm, 
Thereafter, this cut was screened for density by letting them settle in a column of 
water. 

2. The modification of the design of the rods which serve as end-faces for the floating 
liquid column. The new design uses plexiglass rods, machined from optically 
polished plexiglass stock which has been inspected to ensure the parallelism of the 
two surfaces. Thi9 design allows for fine adjustment of the concentricity and 
perpendicularity of the rotating end-faces. Typical values of the total indicator 
readings for concentricity and perpendicularity are 0.0006 in. (0.0015 cm) and 
0.00009 in. (0.0002 cm), respectively. These pieces can be easily changed and 
realigned in the event the end-face9 get scratched. 

3. Improvements in the electronic control of the stroboscopic xenon flash used to 
obtain successive exposures of the tracer particles. The xenon flash apparatus 
consisted of a modified General Radio Company Type 1531 -A Strobotac. by means 
of external circuitry, the capacitance of the discharge capacitor in the Strobotac 
was increased from 1.15 to 17.15 farads at a charging potential of 800 volts, thus 
increasing the energy in the discharge to 5.48 joules. A control box attached to the 
Strobotac trigger input jack was made with circuitry such that a series of trigger 
pulses could be sent to the Strobotac. The maximum number of pulses allowed was 
8, and any of the 8 pulses could be selectively omitted. The time interval between 
pulses could be varied from 5 milliseconds to 1 second. This control box also 
activated the camera motor drive. Opening a switch opened the camera shutter, 
and after a set tim* delay (to allow time for the damping of vibrations), the flash 
sequence took place. Then, by closing the switch the shutter was closed and the film 
was advanced. 

4. The design of a stereoptical photographic system and method to measure the 
particle positions. A stereoscopic camera attachment, also known as a stereo 
reflector, was designed to enable a single camera to obtain a stero-pair of images 
simultaneously. The basic configuration adopted for the design of the stereo reflec- 
tor consisted of a four-mirror arrangement very similar to that used in the Wheat- 
stone stereoscope. 101 

5. The design of a position marker to be used in dark-field illumination for the 
alignment of the optical system. A length of glass optical waveguide, with a core 
diameter of 0.0020 cm attached to an x,y,z micropositioner, was used to create a 
point source of light, with known position with respect to the rotatable end-faces. 
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A photograph of the experimental apparatus is 9hown in Figure 1. The details of the design 
and operation of this system and of the data acquisition method are described in Appendix A. 
The analysis used to calculate the coordinates of a particle from its image location on a pair of 
stereophotographs is described in Appendix B. 


76 


Arthur I) Little Inc 



IV. RESULTS 


Flow Survey 

The position (subject to the qualifications stated in Section III) of the surface separating 
the Proudman-Taylor cells associated with the top and bottom counterrotating end-faces was 
observed and recorded for the test conditions listed in Table 1. 

Dimensional analysis was used to correlate and help interpret the results of the survey of the 
characteristics of the flow within floating columns when induced by counterrotating end-faces. 
The combined action of gravity, inertia, viscous and surface tension forces determines the 
velocity fields within the isothermal floating zone systems under study. Combining this knowl- 
edge with the application of dimensional analysis establishes a convenient set of six independent 
dimensionless variables which determine the flow field: 


• Bond Number, 


pg a i ® 


• Reynolds Number, N R 


P aicoj 


• Weber Number, N w = 


pa? co? 


• a 2 

a i 


• A 

a i 


• co i 
CO 2 

where: 

a = surface tension coefficient of the liquid (dyne/cm), 
p = density of the liquid (g/cm 3 ), 
p = dynamic viscosity of the liquid (dyne-sec/cm 3 ), 
g = gravitational constant (cm/sec 2 ), 
co! = angular velocity of lower end-face (rad/sec), 
co 2 = angular velocity of upper end-face (rad/sec), 
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2 = length of liquid column (cm), 

a x - radius of liquid column at lower end-face (cm), and 

a 2 = radius of liquid column at upper end-face (cm) 

Flow behavior for only a limited set of the dimensionless moduli were investigated. The flow 
characteristics within a short cylinder (2/ai = 1 and = 1) under conditions of counterrota- 
tion were examined in more detail than others. The short cylinder was chosen to reduce the 
relative distortion of the shape of the column due to the effects of gravity and, in addition, to 
provide an approximate low Bond number simulation of a lower gravity field acting on a larger 
prototype system. It follows that the position of the separating plane under the limiting condi- 
tions specified can be given by 


A. = fY^i N N N ) 
h co 2 ’ B ’ w ’ R 

where 

A = distance of the flow-separating plane from the geometric midplane of the cylin- 
drical column. By convention, A is positive for the flow-separating plane above the 
midplane and negative when it is below the midplane; 

h = 2/2 = half length (or half height) of the column. 

Results of selected tests which are characteristic of ail follow. 

Figure 2 shows the results of tests 20 through 24 on the displacement of the flow-separating 
surface from the midplane of a water column having a radius of 0.127 cm and a half-height-to- 
radius ratio of 0.5, with the rotary speed of one end-face Fixed at 1, 2, 5, 10, and 20 rev/sec. A/h is 
plotted versus the ratio of the fixed speed of one end-face to the varying speed of the other, which 
is greater and in counterrotation. The top half of the figure illustrates data obtained with the top 
end-face speed fixed; the bottom half illustrates results with the bottom end-face speed fixed. In 
either case, the plane of separation moved toward the end-face having the smaller rotary speed as 
the speed of one end-face was increased over the other. As the velocity ratio approached zero, 

| A/li | approached unity for low Bond number cases. Also the speed ratio corresponding to any 
selected non-zero value of A/h approached unity as the fixed speed increased. Because the test 
results shown in Figure 2 were all obtained at a relatively low Bond number, the displacements of 
the fiow-separating surface were nearly symmetrical about th* abscissa. 

The uncertainty in the measurement of A/h was estimated to be ±0.05, based on the 
smallest division of the rule in the microscope objective used in taking these measurements. 
Figure 3 shows results from two independent determinations of test number 16. The maximum 
difference between replicate point was 0.15; the average difference was 0.07. 

The effect of increasing Bond number on the displacement of the flow separating surface is 
illustrated in Figure 4. As the Bond number was increased, the static column increasingly 
distorted in the form of a Grecian urn, the position of the flow-separating surface at equal 
counterrotation was displaced upwards, and the pattern of A/h vs. angular velocity ratio became 
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increasingly unsymmetric. Also, at the same Bond number, the position of the flow-separating 
surface shifted monotomically upward with increasing equal counterrotation speed. For instance, 
in tests 16 through 19, A/h at cj 1 /cj 2 =1.0 increased from 0.1 to 0.4. 

Figures 5 and 6 illustrate the displacement of the flow-separating surface with the decreas- 
ing angular velocity ratio that took place with an increase in Weber number at the same Reynolds 
number and with an increase in Reynolds number at the same Weber number, respectively. 

For test number 8 on Figure 5 the surface of separation between Proudman-Taylor cells was 
uncharacteristically below the midplane. It was found that the rods used for this test had slightly 
different radii. The larger radius end-face was on top; the actual radius ratio was a-i/a, = 1.02. 
The larger radius at the top resulted in a greater pumping action of the top end-face in respect to 
the bottom and this, coupled with the change in the shape of the cylinder, drove the position of 
the surface of separation between cells toward the bcHom end-face. 

Figure 7 shows the progression of the position of the flow-separating surface with unequal 
counterrotation up to the onset of the first C-mode instability. The first C-mode occurred at the 
condition correlated by Equation (2) given in Part I, where the rotary velocity limit was that 
reached by either the top or bottom end-face. 

In Figures 8 and 9, values of the arithmetic average of the angular velocity ratio for | A/h | = 
0.65 (this average value will be referred to as V 66 ), as measured in all tests in which h/a was 0.5, is 
plotted against the Reynolds and Weber numbers (based on the fixed and smaller end-face rotary 
speed) pertaining to each test. Figure 8 shows the values of all tests (1-24, except 1, 8, 1 1, 17, 18, 
and 19) for which the Reynolds number is less than 200 and the Weber number is less than 1.0. In 
Figure 9, the scale is changed to a maximum Reynolds number of 800 and maximum Weber 
number of 4.0 to show results for tests 1, 8, 11, 17, 18, and 19. 

For the limiting set of circumstances, shown in Figures 8 and 9, the dimensionless functional 
relationship that applies reduces to: 

V ee. angular velocity ratio = f (N R , N w , N B ). 

In these plots, the average values V are used to reduce the influence of gravity on the 
result. Otherwise the separate influence of Bond number, which ranged between 0.109 and 0.436, 
have been neglected in this display. In any case the maximum difference between the values of 
the velocity ratio of corresponding to A/h = 0.65 and A/h = -0.65 was 9%, with the average 
difference being 4 c c . Accordingly, the assumption that the angular velocity ratio at a fixed value 
of h/a was a function of the Weber and Reynolds numbers only, which is inherent to the plot, is 
reasonable for the relatively low Bond number range tested. 

Further Observations of Surface Swirls 

F * a column of water with end-face radii und height equal to 0.254 cm., the critical speed at 
which swirls first appeared under conditions of equal counterrotation was found to lie between 
3.70 and 3.85 rev/sec. This range of sj ds corresponds to a Reynolds number (based on the end- 
face radii) range of 150 to 156. The pai icular value in this range was found to be influenced by the 
shape of the column as dependent on the variations in zone volume that resulted from filling and 
evaporation. The lower value corresponded to the more obese columns. 
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Both three swirls at a 120-deg separation or two at a 180-deg separation were observed in 
this column for tests at Reynolds numbers ranging between 150 and 4C0, three swirls being the 
more stable state. For tests on a longer column (h/a = 0.75) made at Reynolds numbers between 
150 and 400, the more stable state was found in the two swirls. At Reynolds numbers greater than 
approximately 400, the exact number of swirls could not be distinguished; at these high rotational 
speeds, the flows became visually confusing. Observations in which strobe illumination was used 
did not yield additional information. 

For tests on a water column having end-face radii and a height of 0.127 cm, the onset of 
swirls occurred at equal counterrotation between 14.3 and 16.7 rev/sec. These speeds correspond 
to a Reynolds number range of 145 to 164. For this smaller zone, both the onset of swirls and the 
transition from the state of two to three swirls were more sensitively dependent on slight 
variations of the column shape than for the larger zone. Tests on this column and one having a h/a 
ratio of 0.75 gave similar results as regards the preferred number of swirls. 

The centers of the swirls were observed to precess around the circumference of the column. 
The speed of precession was a function of the angular velocity ratio of the two end-faces. At a 
velocity ratio of unity, the center of the swirls do not precess, but remain fixed in space. As the 
speed of one of the end-faces was increased, keeping the other constant, the swirls moved toward 
the end-face having the smaller rotating speed; moreover, the precessing speed increased and was 
in the irme direction as the rotation of the lower speed end-face. At a sufficiently high-speed 
ratio, the swirls collapsed as they encountered the Ekman layer next to the slower moving end- 
face. 


Velocity Measurements 

Velocity measurements for a column under conditions of isorotation were made in order to 
test the accuracy and precision of the technique on a well defined and understandable circum- 
stance. Under these conditions, the fluid column rotates as a solid body; hence the particle 
pathlines are circles on a plane perpendicular to the axis of rotation and have tangential velocities 
proportional to radius. 

Figure 10 is one of a pair of stereoscopic photographs of a water column of radius - 0.254 cm, 
height = 0.254 cm, under isorotation at 2.0 rev/sec. A six-flash sequence was used, with time 
intervals between successive flashes of .005, .005, .005, .015, and .005 second. Particle images used 
to obtain quantitative position and velocity data are identified. As discussed in Appendix A, 
direct visual interpretation of these photographs may be misleading for it is complicate' 1 by the 
effects of refraction at the liquid-air interface. 

Figure 11 is one of a pair of stereoscopic photographs of the column shown in Figure 10, in 
which an incandescent lamp was used as a light source for an extended exposure (2 sec) From 
this photograph, the radius of the liquid column was determined as a function of the height, and 
these data were used to obtain a fourth-order polynomial expression of the column shape by a 
least squares fit executed by a computer program (as discussed in Appendix B). For this 
particular photograph, a typical case, the shape was described by the equation: 

r = -1.064 Z 4 f 4.142 Z 3 - .1355 Z 2 -.0667 Z 4- .2576. 

The standard deviation of the fit io the experimental data was 5 /j m. 
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The positions of the particles identified in Figure 10 were calculated i sing a computer 
program based on the analysis ’ cribed in Apper lix B. Table 3 shows the calculated values of 
X, Y, Z, r, Mnd co. Figure 12 is a plot of the particle positions in the X-Y plane. In this figure, arc* 
of radius equal to the mean value of r for all flash exposures of each particle are shown, as well as 
the standard deviation about the mean values. The standard deviation for these si> particles 
ranged between 7 and 26 the average value being 11 ^m. The uncertainty interval fur X, Y, 
Z, r, and 0, calculated by error analysis (see Appendix B), are shown in Table 4. The average 
value of Ar calculated fur these particles is ±17 miu. Hence, the deviation of the particle positions 
from a circle is due to the propagation of the random error in the measurement of the particle 
positions on the negatives. 

Note the proximity between particles number 1 and 2 on Figure 12. The average of the 
distance between these two particles for the four flash exposures was 54 ±2^m. The particles used 
in this experiment ranged in diameter from 30 to 45 ^m. These two particles must then be 
clumped together. Computing the vector between these two particles indicated that the two 
particles were traveling as a dumbbell with its axis oriented almost parallel to the Z axis. 5 his 
case serves to illustrate that the resolution of the photographs is better than 54 

The mean of 19 measured values of to is 1.95 rev/sec with a standard deviation of 0.14, n ^his 
mean value is within 2.5 ( c of the known value of 2.00 rev/sec. 

Velocity measurements for a water column of radius = 0.254 cm, height = 0.254 cm, under 
equal counterrotation of 2.5 rev/sec were also made. This circumstance corresponds to test 14 (see 
Table 1), For this test, ilK rotary speed of the end-faces was below that necessary for swirls to 
appear and, therefore, an axisymmetric Flow field was expected. A six-flash sequence was used, 
with time intervals between successive flashes of 0.020, 0.020, 0.020, 0.060 and 0.020 second. The 
calculated values of position and velocity are shown in Table 5. Figure 13 is a plot of the particle 
positions in the r-Z plane and in the X-Y plane. 

The data presented in Figure 13 indicate that the flow follows a simple Proudman-Tavlor 
cellular pattern. The particle path lines wrap on surfaces having the approximate form of toroids 
associated with the ring vortices induced by the top and bottom counterrotating end-faces. The 
dotted line shown is the position of the surface of separation between cells. Although the rumber 
of data points in the top cell is very limited, they are consistent with the ceil pattern described 
and the cell separation plane shown. 

Motion Picture Film 


High-speed moving pictures were taken ot th t circulation within rotating, floating zones as 
revealed by the motion of tracer particles. A steady high-intensity halogen lamp in a dark-field 
arrangement and a 16 mm Ariflex camera operating at 45 frames per second were used for the 
filming. A bellows focusing attachment and a 105 mm focal length Micro-Nikon lens provided a 
magnificatior of 1.0. Kodak Tri-X reversal film was used. 
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The edited film runs for 20 minutes and shows a water column of 0.254-cm radius and 0.254- 
cm reight under conditions of isoiotation, single-end rotation, and counterrotation. Highlights of 
the film include: 

• Part I first shows the column under isorotaticn at 2, 5, and 10 rev/sec. 

• Then, starting at an isoiotational rate of 20 rev/sec, the speed of isorotation was 
increased up to failure in the first C-mode at 43 rev/sec. 

• Part II shows the circulation patterns with top-end rotation of 2 and 10 rev/sec, 
followed by a sequence starting with the column at rest and gradually increasing 
the speed of the top-end through the first C-mode up to failure at the 2nd C-mode. 

• Part III shows results for equal counterrotation at 2, 5, and 10 rev/sec. At 5 and 10 
rev/sec surface swirls are present. The next segment shows the column at equal 
counterrotation starting at 20 rev/sec and thereafter with gradually increasing 
speeds up to the maximum speed of the motors (45 rev/sec). Under thse condi- 
tions the column does not fail. 

The last four segments of the film show the circulation patterns with unequal counterrota- 
tion: 

• First, with the bottom-end rotating at 2 rev/sec and the top-end speed gradually 
increasing from 2 to 20 rev/sec, the., decreasing to 2 rev/sec. 

• A second segment shows the counterpart with the top-end fixed at 2 re* /sec and the 
bottom end varying from 2 to 20 rev/sec and then back to 2 rev/sec. These 
conditions "early correspond to these of test 13 of the flow survey (see Table 1). 

$ The conditions sbotvn on the last two segments of the film correspond exactly to 
those of test ' 5 o. Tie flow 3 urvev. For these conditions swirls are present and the 
film illustrates the characteristics of the skirls already discussed. 


8 : 


i 


Arthur [)!jtt!e Inc 



V. CONCLUSIONS AND RECOMMENDATIONS 


A method for making quantitative measurements of the velocity distributions within small, 
floating columns of liquid using a low Bond number simulation of low-gravity conditions and 
optical-tracer particles has been developed to tne status of proven utility. A survey of the flow 
fields in cylindrical, isothermal liquid columns induced by equal and unequal counterrotation of 
its end faces has been made to provide a background of information for the further experimental 
and analytical investigations needed to make the float zone process frilly understandable and 
predictable. The results of the flow survey revealed the general nature of flows that must be 
predicted by developing analytical constructions and, of particular interest, they showed the 
existence of three-dimensional flows in an axisymmetric geometry and driving condition. 

The velocity measurement technique should now be applied more widely to the purposes for 
which it was developed. Until now, the theoretical predictions of the behavior of floating fluid 
zones and the experiments on them have not been tied together in a rigorous and unambiguous 
manner, due, in part, to ths limitations of both. This situation should — and can — be rectified 
in a cooperative program of experimentation and theoretical analysis. The velocity measurement 
technique should be utilized to provide a data base against which the predictive methods being 
developed with NASA’s support can be tested and improved in an interactive program. 

Experimentation should be extended to the quantitative examination of flow fields within 
heated zones. Moreover, the present capabilities should be expanded by the development of a 
method for extracting the most amount of information on the flow field from photographic records 
of a limited number of particle streaks. The method for extrapolation and interpolation of 
velocity data on a few particles should incorporate the physical principles of conservation of mass 
and angular momentum that prevail and be implemented by computer calculation. Finally, 
minor modifications to the configuration and method of mounting the parts of the optical system 
should be made in order to make the measurement system e&siei to use in extensive application. 
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TABLE 1 


TEST CONDITIONS FOR THE FLOW SURVEY 




a 

h 

u> fixed 

_Nb 

N w 

n r 

Test No. 

Fluid 

(cm) 

a 

(rev/sec) 


(based on co fixed) 

1 

Silicone Oil (3 C.S.) 

.127 

.50 

3 

.369 

.034 

10.1 

2 

Silicone Oil (3 C.S.) 

.127 

.50 

6 

.369 

.136 

20.3 

3 

Silicone Oil (3 C.S.) 

.127 

.50 

10.8 

.369 

.444 

36.6 

4 

Silicone Oil (3 C.S.) 

.127 

.50 

15 

.369 

.850 

50.7 

5 

Silicone Oil (3 C.S.) 

.127 

.50 

30 

.369 

3.40 

10\0 

6 

Silicone Oil (5 C.S.) 

.127 

.50 

5 

.369 

.094 

10.1 

7 

Silicone Oil (5 C.S.) 

.127 

.50 

10 

.369 

.378 

20.3 

8 

Silicone Oil (5 C.S.) 

.127 

.50 

25 

.369 

2.36 

50.7 

9 

Silicone Oil (10 C.S.) 

.127 

.50 

5 

.370 

.094 

5.07 

10 

Silicone Oil (10 C.S.) 

.127 

.50 

10 

.370 

.378 

10.1 

11 

Silicone Oil (10 C.S.) 

.127 

.50 

20 

.370 

1.51 

20.3 

12 

Water 

.254 

.50 

1.25 

.436 

.0111 

50.7 

13 

Water 

.254 

.50 

1.77 

.436 

.0278 

71.7 

14 

Water 

.254 

.50 

2.50 

.436 

.0555 

101.0 

15 

Water 

.254 

.50 

3.54 

.436 

.111 

143.0 

16 

Water 

.254 

.50 

5.00 

.436 

.222 

203.0 

17 

Wattr 

.254 

.50 

7.07 

.436 

.444 

287.0 

18 

Water 

.254 

.50 

10.00 

.436 

.888 

405.0 

19 

Water 

.254 

.50 

20.00 

.436 

3.55 

810.0 

20 

Water 

.127 

.50 

1 

.109 

.00111 

10.1 

21 

Water 

.127 

.50 

2 

.109 

.00444 

20.3 

22 

Water 

.127 

.50 

5 

.109 

.0278 

50.7 

9? 

Water 

.127 

.50 

10 

.109 

.111 

lOVO 

24 

Wa*er 

.127 

.50 

20 

.109 

.444 

2G3.0 

25 

Water 

.127 

.75 

1 

.163 

.00111 

10.1 

26 

Water 

.127 

.75 

2 

.163 

.00444 

20.3 

27 

Water 

.127 

.75 

5 

.163 

.0278 

50.7 

28 

Water 

.127 

.75 

10 

.163 

.111 

101.0 

29 

Water 

.127 

.75 

20 

.163 

.444 

20.3 

30 

Water 

.127 

1.00 

1 

.218 

.00111 

10.1 

31 

Water 

.127 

1.00 

2 

.218 

.00444 

20.3 

32 

Water 

.127 

1.00 

5 

.218 

.0278 

50.7 

33 

Water 

.127 

1.00 

10 

.218 

.111 

101.0 

34 

Water 

.127 

1.00 

20 

.218 

.444 

203.0 

35 

Water 

.254 

.“'j 

2 

.653 

.0356 

81.0 
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TABLE 2 


PROPERTiES OF FLUIDS USED* 


Viscosity 


Fluid 

Density 

Kinematic 

Dynamic 

Surface Tension 


(g/cc) 

(c-s.) 

(dyne-sec/cm 2 ) 

(dyne/cm) 

Distilled Water 

1.00 

1 

0.010 

72.6 

Silicone Fluid 

0.90 

3 

0.027 

19.2 


0.92 

5 

0.046 

19.7 


0.94 

10 

0.094 

20.1 


*These are nominal values at 20° C. The assumption that the liquid column is isothermal at this 
temperature is reasonable, for the accuracy to which calculated values of Ng, and Ng are reported 
in Table 1. 


TABLE 3 

CALCULATED VALUES OF PARTICLE POSITIONS AND VELOCITY 
FOR ISOROTATION AT 2.0 REV/SEC 

Particle No. 


Flash Sequence 

X (cm) 

Y (cm) 

2 (cm) 

r (cm) 

6 Degrees 

co (rev /sec) 

1-1 

-.0424 

.1159 

-.0082 

.1234 

110.11 

_ 

1-2 

-.0348 

.1171 

-.0082 

.1222 

106.56 

1.97 

1-3 

-.0277 

.1200 

-.0081 

.1231 

103.01 

1.97 

1-4 

-.0203 

.1224 

-.0081 

.1241 

99.39 

2.01 

2-1 

-.0406 

.1153 

-.0133 

.1223 

109.42 


2-2 

-.0334 

.1165 

-.0133 

.1212 

106.00 

1.90 

2-3 

-.0284 

.1178 

-.0133 

.1207 

102.64 

1.87 

2-4 

-.0192 

.1209 

-.0132 

.1224 

99.03 

2.00 

3-1 

-.0123 

-.1410 

.0404 

.1416 

264.99 


3-2 

-.0203 

-.1383 

.0404 

.1398 

261.66 

1.85 

3-3 

-.0286 

-.1378 

.0404 

.1407 

258.28 

1.88 

3-4 

-.0364 

-.1359 

.0404 

.1407 

255.01 

1.82 

4-1 

.0561 

-.1130 

.0860 

.1261 

296.39 



4-2 

.0491 

-.1196 

.0860 

.1293 

29^.32 

2.26 

4-3 

.0419 

-.1210 

.0860 

.1281 

2r9.09 

1.79 

4-4 

.0348 

-.1224 

.0860 

.1272 

2t 5.85 

1.80 

4-5 

.0116 

-.1287 

.0361 

.1292 

275.13 

1.98 

5-1 

.1246 

-.1705 

.0441 

.2111 

306.15 



5-2 

.1140 

-.1825 

.0440 

.2152 

301.99 

2.30 

5-3 

.1027 

-.1872 

.0440 

.2135 

298.75 

1.80 

5-4 

.0914 

-.1969 

.0440 

.2171 

294.92 

2.13 

61 

.0464 

.2542 

-.0920 

.2584 

280.34 



6-2 

,0311 

-.2587 

-.0920 

.2606 

276.86 

1.93 

6-3 

.0161 

-.2600 

-.0919 

.2606 

273.56 

1.84 

6-4 

.0003 

-.2614 

-.0919 

,26 1 4 

270.08 

1.93 
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TABLE 4 


RESULTS OF ERROR ANALYSIS OF SELECTED PARTICLES. 
BASED ON A VALUE OF THE MAXIMUM UNCERTAINTY 
IN THE POSITION ON THE NEGATIVES OF ± .00038 CM 



AX 

AY 

AZ 

Ar 

A 0 

Pfcrtide 

(cm) 

(cm) 

(cm) 

(cm) 

(Deg) 

1-1 

.0002 

.0009 

.0003 

.0008 

.19 

3—1 

.0003 

.0017 

.0003 

.0017 

.13 

4-4 

.0003 

.0017 

.0003 

,0016 

.22 

5-4 

.0004 

.0021 

.0003 

.0020 

.20 

6-1 

.0004 

.0022 

.0003 

.0022 

.10 

Mean 

.0003 

.0017 

.0003 

.0017 

.17 
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CALCULATED VALUES OF PARTICLE POSITIONS AND VELOCITY FOR EQUAL COUNTERROTATION AT 2.5 REV/SEC 


nm 

o’ o d d o 


m 

u c o 


ddo 


8 $i 

odd 


s 

o 


? 88 : 


s = 


!!§ 


382 

PM IN CN 


i 


as §S 8 

o o o o o 


« 8 S t 

555 


5 5 «? 


2 i 

9 


8 §i 

9 9 9 


§' §> §sl§5 


o o o o o 


ss 

CM CM 

9 9 


• 

9 9 9 


§88 


ss n® , 
& 55555 


i§ > §§§ 


ass , 


PM | 
m 1 


ID (O 00 ID Pm 

- - N ffl 

<C V D CN «- 


a m i 

a a 1 

o o 


is? 


I 


a a a oo a 

88 553 

9 9 9 9 9 


£ ' 


5 £ 


0(0 0 
o rv r«. 
o o o 
pop 


O irt CN 

Sss 

9 9 9 


2 88 
*7 a a 

9 9 9 


5 8 £ 8 g 

O *— CN •— O 

o o o o o 


s s 

o o 
o o 


S 8 5 


88 E 

CO o— «- 

o o o 


1 

I 


S 5 (M Oi O 

O ui iS 

nj d p CN 


m in a 

u> — a 

o cn d 


? £ 8 

CN CN (N 


S S 5! 

(N CN CN 


a p*- m cd «■ 

N N IN (\ d 

tfi ift if! ID Mf 


m a 
in in 


S m ' 

cn 


cn o in 
O cn o 


1 


nsn sss 


o to »— op in cn 
^ f o o> o 

CN CN CN CN *- *— 


s; 


? s s 


n - id c 

£ 3 8 8 


£ S 

a o 


a a «- a 

n ^ 

in S In S 

in m re cn 


. a 

I' 03 


3 8 S £ ? S 

CN «r O (O 2 O 
a id a a a O 

CN CN CN CN CN cn 


5 8 S a 8 < 


CN r- Q O 
00 CT) O •— 

CN CN cn P3 


a — in a 
— cn « m 
n n i o n 


r- f- cn m 3 m 
^ q o o cn in 

S3 3 3 3 3 

o © o o o o 


cn cn rt 
op cn cn 
cn mt m 


cn cn oo cn 
5 cn o 2 

3 5 a 

o o o o 
o o d o 


CN (N D *• 
in ^ cd in 
00 05 h N 

o o o o 

doc ; 


Mf cn 

iN CD 
CN •- 


— m 
pi cn 
cn r j 


*N © O © 

co oo cn in 
ill M v a 

o o o o 


5 8 


OI ID to lO ID M 

ooainoa — 

CN CN CN CN CN CN 

O O O O O O 


co cn a 
— cn n 
cn cn cn 


cn — *» 
<- m O) m 

CN CN PM CN 

o o o o 


mt -or a 
a> 5 a 8 

o d o o 


*C ^ cd r- uS in rs cD 

I 8 8 8 8 S S 8 

o o d d o o ’’ 


m oo cn in 

8000 

pood 


IRS 


o o o o 
9 9 9 9 


o o 
o o 


88 

o o 


? § rl OI 
8 8 8 8 
o o o o 


8 8 

o o 


88 

□ o 


in oo 03 a a 

<D CN — 00 Pi 

CNCnCNCNt—t— 

o c o o o o 

ope o'do 


5 8 8 

poo 


cn <D a ID 
cn in m — 

8 8 8 ° 

9 9 9 ° 


Mt *- m rv 
■j a o a 

8 8 8 8 

p o o o 


E 

> 


O a *- O n- q 
cp cn v a — o 
PM cn cn — ~r *- 
o o o o o o 
9 q pood 


oo a a 

288 


o o o o 
o d o o 


a a »- a 
Pi PM 00 Mf 
CN CN — — 

o o o o 
d o o o 


88 

o o 


— in id o 

CD CN m n 
9 cr rt n 


3 S 

o o 


5 3 

o o 


o o o o o o 


6 o o 

odd 


S CO ffl s 

in o ig 
9 9 9 9 


5 g i § 

9 9 9 9 


E 

u 

X 


*• n o oo oo m 

pi a S 8 5 ? 

o o o o o o 

o p p pop 


8 2 5 

9 9 9 


p o«- n 

v 5 R n 
© o o o 
9 9 9 9 


# <“ CN N 
O 90 IO PI 
S N N N 

o o o o 
d o o o 


8 a 
8 


d o o o 
o o o o 


00 5 
a o 

pi pi 


cn o 3 oo cn n 
a a a a a a 
in c o <n o pm 
o c o o — 
o o o o o o 


Id in 3 8 

S 3 8 8 

o o o o 


S3! 


»- cn a ^ in p 


— cn a 
fit CM N 


- CN Cl M 

rt pi cn rt 


•- a 
a a 


•— a a ip cc 
d o d d o o 


a a a a 




88 


» 


Arthur D I jttlc Inc 



SPINDLE SPEED 
CONTROL UNIT; 



EXPERIMENTAL SET-UP FOR INVESTIGATION OF FLOW WITHIN FLOATING ZONES 










Angular Velocity Ratio 

FIGURE 2 EFFECT OF FIXED ANGULAR VELOCITY ON A 
AS A FUNCTION OF ANGULAR VELOCITY h 
RATIO 


90 


Arthur D Little; Inc 


i 









=r|> 


1.0 



Angular Velocity Ratio 

FIGURE 4 EFFECT OF BOND NUMBER ON-j^AS A FUNCTION 
OF ANGULAR VELOCITY RATIO 
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FIGURE 5 


Angulsr Velocity Ratio 

EFFECT OF WEBER NUMBER ON^- AS A FUNCTION OF ANGULAR 
VELOCITY RATIO, AT A CONSTANT RFYNOLDS NUMBER 
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FIGURE 6 EFFECT OF REYNOLDS NUMBER ON^-AS A FUNCTION OF ANGULAR 
VELOCITY RATIO, AT A CONSTANT WEBER NUMBER 


94 


Arthur H Little Inc 






.1 .2 .3 .4 .5 .6 .7 .8 .9 1.0 

Angular Velocity Ratio 

FIGURE 7 ^ AS A FUNCTION OF ANGULAR VELOCITY RATIO FOR CONDITIONS 

APPROACHING THE 1ST C-MODE INSTABILITY 
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FIGURE 9 PLOT OF V .65 AS A FUNCTION OF REYNOLDS AND WEBER 
NUMBERS FOR TESTS 4, 5, 8, 11, 17, 18, 19, 20, 22, 23, AND 24 




FIGURE 10 ONE OF A STEREOPAIR OF PHOTOGRAPHS FOR 
ISOROTATION AT 2.0 REV/SEC 



FIGURE 11 PHOTOGRAPH USED TO MEASURE THE SHAPE OF 
THE FREE SURFACE OF THE COLUMN 
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APPENDIX A 


DESIGN CONSIDERATIONS FORA 
STEREOPTICAL SYSTEM AND A METHOD FOR 
TAKING VELOCITY MEASUREMENTS 


DESIGN 

The design of the stereoptical system is illustrated in Figure A-L The basic principle 
adopted for the design of the stereo-reflector is a four-mirror arrergement very similar to that 
used in the Wheatstone stereoscope. 111 

The stereo-reflector consists of plane mirrors m,, m„ and m 4 . The two outer mirrors, m, and 
m„ are arranged at an angle less than 45 deg to the axis of the lens. Between those two mirrors, 
the inner pair is positioned so that the reflecting surfaces are at right angles to each other and at a 
45-deg angle to the axis of the lens. The inner vertical edges of m» and m« meet on the optical axis 
and, as a result, each of the two mirrors occupies one half of the field of view. 

A “left-eye” view of the subject is formed by rays, such as L, which after reflection at the 
surface of m, are then reflected in the direction of the lens by m,. The image is brought to focus 
at E . Likewise, a “right-eye” view is formed by rays, such as R, which after reflection by the 
mirrors m a and m« are brought to focus at r. 

The images r and E can be thought of as being formed by two cameras rotated by an angle 
0, and 0, to the left and right of the optical axis, as shown by the dotted lines. The positions of the 
two virtual cameras and the actual camera are related by the following geometrical relationships: 


0, =2(0, -45°) -a, 

(A-l) 

0 2 = 2(0 2 - 45°) -a 2 • 

(A- 2) 


where 

0i = the angle between the optical axis of camera 1 and the y-axis of the object 
coordinate system, 

#i = the angle between the plane of mirror m, and the perpendicular to the optical 
axis of the actual camera, and 

a, = the angle between ray L and the optical axis of the actual camera. 


(1) Applied Optics and Optical Engineering (Klngslako, R. ed.). Academic Prats, Vol. II., IMS, p. 89. 
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A similar nomenclature with subscript 2 applies to ray R: 


Oj = tan' 1 

d, 

FL 

(A-3) 

a 2 = tan' 1 

FL 

(A-4) 


where 


dj = the distance from the center of the image at G to the center of the film plane, 

d, = the distance from the center of the image at r to the center of the film plane, 

FL = the film-to-lens conjugate distance, 

OL = the object-to-lena conjugate distance. 

Assuming the camera lens is thin, then: 


FL = ( 1 +M)/ 



(A-5) 

(A-6) 


where 


M ~ magnification, and 
f = the focal length of the lens. 

The first step in the design of the system is the selection of values of f and M. Both 
theoretical and practical considerations must be taken into account. Some of the theoretical 
aspects of photographing a small particle in dark-field illumination have been treated in detail in 
previous work. {,) For smooth spherical particles illuminated by light coming from one direction, 
assuming that the light can be treated as consisting of plane waves and that it is scattered 
isotropically by the particle, the camera lens will produce a diffraction pattern of the glint of light 
it receives from the particle. If this diffraction pattern is the classical Airy disk, the diameter of 
the central disk is given by: 

D a = 2.44XF (A-7) 


where 


X = the wavelength of the light, and 
F = the f-number or relative aperture of the lens. 

(2) Fowls, A. A., et al M Float Zona Processing In a Weightless Environment, NASA CR-2768, 1976. 
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The F-number or relative aperture of the lens iB defined as: 


F = 


1 

D 


(A-8) 


where 


f is the focal length, and 
D is the diameter of the aperture. 

In macrophotography, the f-number is defined as either: 

F-O+M)/ /(1+M) (A-9) 

D or MD 

whichever is smaller. This formulation represents the so-called bellows correction. 

The diameter of the resolution element, AX, in the object plane is found by dividing the 
diameter of the resolution element at the image by the magnification M. Thus: 

(A - ,0) 

The depth of field, calculated by geometrical optics, is given by: 

AY= 144 X (s) 1 - <A ' n> 

Equations (A- 10) and (A-ll) show that the depth of field and the size of the resolution element in 
the object plane depend only on the ratio F/M (for light of e given wavelength). We would like AX 
to be small (high resolution) and AY to be large (larg* depth of field), but these are conflicting 
requirements. 

In practice, the higher limit of the F-number is determined by the following parameters: 
particle size, particle reflectivity, energy output of the illumination system, and film speed. The 
ease with which the images can be enlarged and accurate measurement of particle position be 
made is an added consideration that must be taken into account in choosing the magnification. 

Taking all the above factors into account resulted in the following design specifications: 

• a 105-mm, f/4 Micro-Nikkor lens, 

• M = 0.80, 

• F * 24.7, 

• Nickel-coated, hollow-glass spheres 30-45 in diameter, 

• Xenon flash apparatus consisting of a modfied General Radio Company, type 
1531-A Strobomer, with discharge energy of 5.48 joules, and 

• Kodak Plus-X 125 ASA film. 
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For these conditions, adequate contrast, a depth of field of 0.13 cm (on the order of the 
radius of the liquid column), a resolution of 42 nm (of the order of the particle diameter), and an 
image size of 0.2 by 0.4 cm which could be enlarged to useful size with a microscope at 20X were 
achieved. 

The design specifications of the stereo reflector were tailored to the specific value of f-105 
mm and M=0.80. A schematic diagram of the device is shown in Figure A-2. For mirrors mi and 
m„ 25.0 mm square, first-surface mirrors (Melles-Griot #02 MPG 001) were used. A right-angle 
prism (Melles-Griot #01 ORS 020/011) with silvered legs was used for mirrors m, and m 4 . 
Machined parts were used to hold the two mirrors and the prism in place. Four sets of these parts 
were made; they are interchangeable and provide angles B x and 0 3 of 51, 53, 55, and 57 deg. These 
correspond to convergence half angles of 8.0, 12.0, 16.0, and 20.0 deg., respectively, for f= 105 mm 
and M = 0.80. 

The following procedure was followed in aligning the system for making quantitative 
measurements: 

Step 1 — The alignment of the two plexiglass rods is verified, with minor adjustments 
made if necessary; 

Step 2 — The X, Y, Z, micropositioner is calibrated such that the position of the end 
of the optical fiber is known with respect to the center of the end-faces; 

Step 3 — The magnification is set by adjusting the bellows attachment and measured 
from photographs of a scale; 

Step 4 — The camera is aligned on the optical bench, such that the optical axis passes 
through the center of the column and is orthogonal to the coordinate system 
of the micropositioner; and 

Step 5 — The stereoreflector is aligned such that both stereo images are in focus and 
d, - d s . 

The system can then be used to obtain stereoscopic photographs of the liquid column from which 
quantitative measurements can be made. 

IMAGE INTERPRETATION 

In interpreting the image of particles within the liquid column, refraction effects at the 
water-air interface must be taken into account. Figure A-3 shows the effects of refraction at the 
front surface of a water column. The positions of particles within the column are significantly 
distorted by refraction; only particles very close to the font surface of the column are not 
affected. Particles in the zone labeled “shadow zone” in Figure A-3 will be invisible to the 
observer, and worse still, the very existence of this shadow zone will not manifest itself to the 
observer at all. Rays emitted by particles in the zone labeled “double image zone” will emerge 
near the edge of the cylinder, giving the impression that the contents of the entire cylinder are in 
view. For particles in this zone, we see a double image. 

Figure A-4 shows the overlap between the shadow zones and double image zones for the two 
cameras at a convergence half angle of 12 deg, focussed about the center of a liquid column of 
radius = 0.127 cm. Also shown are the areas of best focus for each camera for a depth of field of 
±AY = 0.13 cm. The result is 14 distinct zones in the liquid column that must be taken into 
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account in interpreting a stereoscopic pair of particle .u~ges within the water column. ThiB 
makes the visual interpretation of the stereo photographs a rather difficult task. The areas for 
which useful information can be obtained within the liquid column are 1, 2, 3, and 4 as labeled in 
Figure A-4. 

DATA EXTRACTION 

An Opto-Metric Tools Model 1041 toolmaker’s microscope was used to make measurements 
directly on the 35 mm negatives. The negative is mounted on a glass slide which is then aligned, 
such that the X ane Z directions on the negative are parallel to the X and Y directions of the 
microscope as shown schematically in Figure A-5. 

Since the camera is focussed about the center of the column, and the optical axis passes 
through the center, the outline of the column will exhibit no parallax (assuming the shape of the 
column is axi-symmetric. Therefore, the edge can be used to determine the radius as a function of 
the height. These data are used in correcting the calculations of particle positions for the effects of 
refraction. This is discussed in Appendix B. 


Measurement of the position of the particle images is made as follows. Referring to Figure 
A-5, one sees that the origin of the film coordinate system is determined by: 


X = la : b ' + a 
0 2 

(A- 12) 

31,(1 l Al 

Z a 1 C_ d ! +d . 
0 2 

(A-13) 

The coordinates of point e are determined by: 

X e = X 0 -e x 

(A- 14) 

and 

Z =e -Z . 

e x o 

(A- 15) 


where 


a, b, c, d, e„ e, are positions expressed in the microscope table coordinate system. 

The reproducibility in the determination of the location of a single point using the toolma- 
ker’s microscope is ±.00005 in. (±1.3 jim). The maximum uncertainty in the measured value of 
X. and Z, is three times the value of the uncertainty in the position of a single point or ±3.8 ^m. 
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FIGURE A-1 SCHEMATIC DIAGRAM OF STEREOPTICAL PHOTOGRAPHIC SYSTEM 
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FIGURE A-4 CONSTRUCTION SHOWING ZONES OF DIFFERENT PARTICLE IMAGE 
QUALITY WITHIN A CYLINDRICAL COLUMN OF WATER 
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APPENDIX B 


CALCULATION OF PARTICLE POSITION AND VELOCITY 


METHOD 

The position of a particle in the liquid column is determined from its location on a 
pair of stereoscopic photographic images. The analysis is based on ray traces in the optical 
system shown in Figure B-l . For each stereoscopic pair, a ray is traced starting at the image 
position, passing through the center of the lens. The points of intersection of each ray with 
the surface of the liquid column is then found, and by applying Snell’s Law, the directions 
of the refracted rays are obtained. The points of closest approach between the two refracted 
rays are determined and their midpoint taken to be the position of the particle. 

Three co-ordinate systems are used in the analysis: the object coordinate system and 
two film plane coordinate systems (f, and f 2 ), as shown in Figure B-l. 

A ray of arbitray lenth originating at a point on film planes f, and f 2 can be written 
as: 


-** -*■ A 

Pi =b 0 i +Pis, 

--k ^ 

(B-l) 

P2 = ^02 + P2 S 2 

where 

— V. 

b 0 ! = base of ray expressed in object coordinates, 

Pi = direction of the ray (unit vector), and 
S! = length of the ray. 

The nomenclature follows similarly for subscript 2. 

(B-2) 

The position of a point on the film plane is transformed from film plane to object CO' 
ordinates by the following rotaion and translation equations: 

bo 1 = Ri b f i + Tj 

— k. 

(B-3) 

bo 2 = ^2 b f2 + 

(B-4) 
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* here 


b 0 1 = base of ray in object coordinates, 
b^, = base of ray in film plane coordinates, 

R, = rotation maxtrix from the film 1 coordinate system to the object 
coordinate system, and 

T, = origin of the film coordinate system expressed in the object coordinate system. 
The nomenclature follows similarly for subscript 2. 

By substituting the coordinate values of bf , , Rj , and T , into Equation (B-3), it can 
be rewritten as: 


( cos0j-sin0, °\/ x »\ /-(FL + OL)sin 0 ,\ 

sin <t>i cos 0i 0 if 0 I + I — (FL + OL) cos 0, J 

0 0 l/Vz,/ \0 / 


(B-5) 


Similarly substituting the values of b f 2 , R 2 , and T 2 into Equation (B-4), it can be rewritten 

as: 


^ ' cos0 2 sin 0 2 0\/x 2 \ / (FL + OL)sin 0 2 \ 

b 02 = sin 0 2 cos 0 2 0^0 J + (FL + OL) cos 4> 2 J 


(Carrying out the matrix multiplication and addition yields: 


(B- 6 ) 


&oi ~ 


boa = 


( Xi C'js0j — (FL + OL)sin 0A 
y, sin 0 , - (FL + OL)cos0, I 

( x 2 cos 0 2 + (FL + OL) sin 0 2 \ 
- x 2 sin 0 2 - (FL + OL) cos 0 2 J 


(B-7) 


(B- 8 ) 


Rays L, and L } , originating at the vectors b 0 1 and b 0 2 , must pass through the center of 
the thin lens. These rays in the object coordinate systems are: 


/— OLsin 0 ,\ 
-i = l — OL cos 0 i j 


(B-9) 
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(B-10) 



The unit direction vectors for rays P, and P 2 are: 


*»- 


Li — bpi 

IL, -tj,l 


^2 * 



(B-ll) 


(B-12) 


The next step is to find the point of intersection of rays P t and P 2 with the surface of the 
liquid column. 

Due to the effect of gravity, the shape of the column is not quite cylindrical. The 
radius of the distorted column is a function of the height (z) and is described by the 
following equation: 

r = f(z) = a n z n + a n l z n ' 1 + . . . + a,z + a 0 (B-l 3) 


where 

r = radius, 

z = height, and 

a’s = regression coefficients of a least squares fit to experimental data. 

The surface of the distorted liquid column is described by the equation: 

g(x,y,z) = x 2 +y 2 - (f(z)] 2 = 0. (B-l 4) 

Substituting the ray Equation (B-l) into Equation (B-l 4) yields the expression: 


g(P,)= [b 0 )x +Pi x sj 2 + j^boiy + Pi y s J 2 — £f(b<n s + Pi 2 S )J 2 (B-l 5) 

Equation (B-l 5) is of the general form g(s) = 0 and can be solved by iteration using Newton’s 
method, which states: 


g (Sj) 

7 <sj) 


(B-l 6) 
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s(s ) 

This is repeated until | . ** n - | < e . 

8< s i+„> 

As g(s) is a continuous single-valued function of X, Y, and Z, its first derivative can be 
expressed by: 


g'(s)= dx + dg^ dy + d£ dz^ 

dx ds dy ds df dz ds 


Differentiating Equation (B-15) yields: 


| = 2 I>* + P.xs] 


dx _ c 
ds 


dy 2 [^ 01 y + Piy S J 


dy /\ 

5T =P ” 


f-™ 


df 


jl = na n z " " 1 +n-la n _ 1 z n - 2 +... + a 1 


n-2 


dz 


dz * 

dT Pl > 


(B-17) 

(B-18) 

(B-19) 

(B-20) 

(B-21) 

(B-22) 

(B-23) 

(B-24) 


To calculate a first approximation of s to be used in Newton’s method, we find where 
the ray intersects a right circular cylinder of radius R whose surface is described by the 
equation: 


x J + y 2 - R J = 0 


(B-25) 
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Substituting the ray Equation (B-l) into Equation (B-25) yields the expression: 

[b 0 ix + PixS] J + tboiy + Pi y s] 2 - R 2 = 0 . (B-26) 

Multiplying out gives: 

l(Pix) 2 + (Piy) 2 Is 2 + {(b 0 ix Pix)+(boi y Piy)ls+t5oix) 2 +(b 0 ,y) 2 ] -R 2 =0 . (B-27) 

Equation (B-27) is quadratic in s and can be solved explicitly. This result is then used as an 
initial value in Newton’s method. 

In practice, the iteration is repeated until lg(s)/g' (s)l < .0001 cm; the final value of 
s calculated was then substituted into the ray Equations (B-l) and (B-2) to obtain the points 
of intersection pu and p 2 j. The normal to the surface of the liquid column is then 
calculated: 



The normal for the surface of Equation (B-l 4) is: 

n = ===== ===. ((x 2 y ! — f (z)f ' (z)) . (B-29) 

y x 2 + y 2 + (f(z)f'(z)] 2 

The ray is refracted as it passes into the liquid. The unit vector q of the refracted ray 
can be determined by applying Snell’s Law, which states that the incident ray p, the normal 
n, and the refracted ray q lie on the same plane, and 

nj sin0, = n 2 sin 0 2 (B-30) 

where 

nj , n 2 are the refractive indices of the two media: 

0j = the angle between p and n, and 

0 2 = the angle between q and n. 

The term q can be expressed in terms of its components relative to n and w, where w is 
the unit tangent vector lying in the plane defined by the incident and refracted rays (plane 
of incidence): 

q = -ncos0 2 + wsm0 2 (B-3 1 ) 
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XJ> 


where 


A A 

a a a , a pxn 
w = n x v and v = -r— - 
sin 8 X 

is the normal to the plane of incidence. 
Hence: 


w = — — -r— |n x ($ x n)] 
sin 0j 


w = 


1 


sin 

1 


w = 


sin 0, 

Substituting (B*33) into (B-31) yields 


[(n*n) p - (j$*n)n] 


[p + n cos 0i ] 


Note that 


A a a sin 02 -■ a a a | 
q = n cos 0 2 + - — — [p + n cos 0j ] 
sin 0 , 

n 2 n 2 


cos 0 j - - n*p 


cos 


s0 2 = ■y/ 1 - sin 0 2 
sin 0, — 1 

n 2 


sin 0 2 = 


n 2 


In x p| 


(B-32) 


(B-33) 


(B-34) 


(B-35) 

(B-36) 

(B-37) 


can then be calculated by combination of Equations (B-34) through (B-37). 

The refracted rays can be expressed as: 

Qi =Pii + qi t, 

Q2 =p 2 i + qj t 2 


(B-38) 

(B-39) 


where t, and t 2 are the length of the rays. 

The next step is to find the value of tf and tf for which the distance D between the two 
rays is a minimum. The distance D is calculated as: 

D= ^(Q, - Q 2 ) • (Qi Oi) (B-40) 
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Substituting Equations (B-38) and (B-39) yields: 


D = }[( pit + q,tj - p 2i - q a t 2 )• (Pn +qiti -p a i -q 2 t a ) (B41) 

Differentiating Equation (B-41) with respect to ti and t a yields: 


3D a r* , a . a , v 

r— = 9i *(Pii + 9i ti - Pai - 9a U ) 
dt, 

(B42) 

3D -a r* , a . a . \ 

= _q 2 .( Pli +qitj - p 2i - q 2 ta) 

Otj 

(B-43) 

Letting the derivatives equal zero gives two equations with two unknowns, tf and tf : 

(<li *9i )tf — (qi *q 2 )t* ~ 9i *(Fii — P 2 i) 

(B-44) 

- (q j *q 2 )t * + (9a *92 )t* = 92 • (Pii - P 2 i) 

(B45) 

Letting: 


(9i-9i) = a, 

(B46) 

— (9l *92 ) = *2 

(&47) 

(92 * 92 ) = 33 

(B48) 

* (Pli — P2i) = hi 

(B-49) 

A -A 

q2 * (p li p2i) “ ^2 > 

(B-50) 

we can write: 


3j tj + ^2 = hj 

(B-51) 

ajtj 3jtj “ bj 

(B-52) 

Solving these two equations gives values for tf and tf : 


* W a3 — bjaj 
a t a 3 - (a 2 ) 2 

(B-53) 

and 

_ b 2 aj — bia 2 
a l a 3 — ( a 2 ) 2 

(B-54) 
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Substituting the values of tf ancHf intoEquations (B-38) and (B-39) yields the value of the 
two points of closest approach Qf and Qf and the position of the particle in the liquid 
column is: 

(X,Y,Z) = y(Qf + Qf) (B-55) 

where Qf and Qf are expressed in terms of their object coordinate system. 

The particle positions are also calculated in cylindrical coordinates: 

r = [x 2 + y 2 ] 1/2 
and 

0 = tan -1 - 
x 


(B-56) 

(B-57) 


ACCURACY 

To calculate the uncertainty in the results X, Y, Z, r, and 0, we express the results as a 
function of the independent input variables Xj , Xj , z, , and z 2 : 


X » f (Xj ,x 2 ) 

(B-58) 

Y = f (x,,x 2 ) 

(B-59) 

Z = f (Zj ,z 2 ) 

(B-60) 

r = f (x l5 x 2 ) 

(B-61) 

e = f (x, ,x 2 ) 

(B-62) 


The uncertainty intervals AX, AY. AZ, Ar and Ad are calculated by the following expressions: 



(B-63) 


(B-64) 


(B-65) 
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(B-66) 



(B-67) 


Equations (B-63) through (B-67) were solved numerically using the computer program to 
calculate the results, while varying each of the four input variables by a differential amount 
equal to 4pm, which amount corresponds to the reproducibility of the measurement of a 
point on a negative using the toolmaker’s microscope. The results of applying this method 
to the measurement of velocities in a case of a column rotating in isorotation is presented 
in Table B-l. 
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TABLE B-1 


RESULTS OF ERROR ANALYSIS OF SELECTED PARTICLES, 
BASED ON A VALUE OF THE MAXIMUM UNCERTAINTY 
IN THE POSITION ON THE NEGATIVES OF ± .00038 CM 



AX 

AY 

AZ 

Ar 

AO 

Particle 

tcm) 

(cm) 

(cm) 

(cm) 

(Ota) 

1-1 

.0002 

.0009 

.0003 

.0008 

.19 

3-1 

.0003 

.0017 

.0003 

.0017 

.13 

4-4 

.0003 

.0017 

.0003 

.0016 

.22 

6-4 

.0004 

.0021 

.0003 

.0020 

.20 

6-1 

.0004 

.0022 

.0003 

.0022 

.10 

Mean 

.0003 

.0017 

.0003 

.0017 

.17 
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